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Abstract 


This  report  develops  a  scoring  or,  equivalently,  validation 
procedure  for  those  multiple  target  correlator- trackers  which,  in 

O 

effect,  form  a  partitionii^  Q  of  accumulated  data  Z  from  multiple 

sensor  sources  into  disjoint  track  sets  and  a  possible  false  alarm 

set.  The  accumulated  data  may  involve  missed  detections,  false 

alarms,  and  may  describe,  typically,  measured  target  positions, 

velocities,  or  various  target  attributes,  such  as  hull  length  or 

even  identification.  One  such  score,  J,  for  a  given  correlator- 

«  0 

tracker  is  defined  to  be  J(Q,Z)  =  -2  log  pr(Q|Z),  a  form  related 

to  the  posterior  distribution  function  of  partitionings  of  data. 

Use  of  J  as  a  measure  can  be  justified  from  both  information  and 

decision  theoretic  viewpoints.  In  particular,  minimal  J  is 

achieved  by  the  posterior  maximum  likelihood  estimator  of  Q. 

0  eo 

^CQ,Z)1Q),  a  measure  of  cross  entropy,  is  minimized  among  constant 
'  0  00  0 

partitionings  by  Q  =  Q.  A  closely  related  score  J'(Q,S)  =  -2  log 

O 

Pt(z|Q)  has  similar  properties  to  J  and  is  computationally  more 
convenient.  The  emphasis  of  this  report  is  on  J'. 

When  a  linear  Gauss-Markov  tracking  and  observation  model  is 
assumed,  and  false  alarms  and  attribute  information  are  modelled  by 
Gaussian  processes,  J'  is  shown  to  be  a  sum  of  relatively  simple 

O  O 

computable  terms.  The  distribution  of  (J'(Q>Z)|Q)  under  the  above 

assumptions  is  that  of  a  chi-square  random  variable  plus  a  constant. 

Thus,  3'  is  a  natural  measure  of  the  goodness  of  fit  of  a  correlator- 

tracker's  output  to  the  data  it  is  operating  on. 

Several  correlator-trackers  can  be  ordered  with  respect  to  overall 

relative  accuracy  through  use  of  their  scores  in  a  weighted  sense  - 

depending  on  the  prior  decision  probabilities  and  the  decision  costs 

involved.  However,  it  is  also  shown  that  type  I  and  type  II  decision 

errors  are  difficult  to  compute,  since  distributions  involving 
0  00  0 

C(J'(Q,Z)  -  J'(Q,Z))jQ)  apparently  cannot  be  obtained  in  simple 
form,  in  general. 
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INTRODUCTION 


^^Jrhis  report  is  the  first  part  of  a  two  part  study  concerning 
performance  measures  and  scoring  procedures  for  the  multiple 
target,  multiple  sensor  correlation  problem.  In  this  document^a 
methodology  is  developed  for  determining  how  well  the  outputs  - 
in  the  form  of  partitionings  of  the  data  -  of  a  given  correlator- 
tracker  match  the  actual  track  sets  and  false  alarms  that  are 
present. the  second  part  of  the  study,  numerical  examples 
will  be  presented  for  simulated  real-world  situations,  including 
ocean  surface  and  stibsurface  scenarios.  These  examples  should 
both  illustrate  the  potential  use  of  the  scoring  procedures 
developed  here,  and  lead  to  various  sensitivity  analyses.  In  the 
latter,  trade-offs  are  expected  to  be  established  between  a  number 
of  quantities  of  interest,  including,  background  shipping  densities; 
types  of  evasive  target  motions  and  group  formations;  number  of 
(true)  targets  of  interest;  false  alarm  rates;  incoming  sensor 
measurement  rates;  averaged  measurement  errors;  and  correlation- 
tracker  design  -  illustrated,  e.g.,  by  the  number  of,  or  average 
length  of,  track  sets  formed. 

The  study  is  part  of  the  on-going  Correlation  Handbook  Project 
at  the  Naval  Research  Laboratory  (Code  7932) ,  under  the  aegis  of 
NAVELEX.  The  results  of  this  report  extend  and  supersede  much  of 
the  earlier  work  presented  in  Ref.  [Ij. 

A  performance  (or  validity)  measure  of  a  given  correlator- 
tracker  determines  -  through  its  numerical  outcomes  -  how  well  the 
outputs  of  that  correlator-tracker  estimate  the  true  target  - 
false  alarm  situation,  through  its  partitioning  of  the  observed 
data  into  track  sets. 
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The  purpose  of  this  report  is  to  present  a  scalar  perfor¬ 
mance  measure  which:  is  applicable  to  a  wide  class  of  correlator- 
trackers;  is  relatively  simple  in  form,  suitable  for  real  world 
operational  implementations  for  use  in  ocean  surface,  subsurface, 
air,  or  ground  surveillance,  or  in  any  combination  of  these 
scenarios;  and  is  justified  by  a  mathematical-logical  basis. 

(For  additional  explanation  of  terms  used  here,  see  the 
following  section.) 

The  correlator-trackers  treated  here  operate  on  data,  re¬ 
ceived  by  possibly  several  sensor  systems,  which  is  in  the  form 
of  reports  that  may,  unknown  to  the  observer,  represent  targets  of 
interest  or  false  alarms.  The  class  of  sensor  systems  considered 
is  general  and  is  essentially  only  limited  by  the  restriction  that 
the  time  and  sensor  source  of  each  report  be  known.  Reports  may 
contain:  (1)  geolocation  target  or  false  alairm  measurements,  such 

as  observed  positions  or  velocities,  which  typically  come  from 
bearing,  range,  and  related  sensors;  and/or  (2)  possible  non¬ 
geolocation  target  attribute  information,  such  as  observations  of 
hull  length,  number  of  sensors  on-board,  emitted  frequency  infor¬ 
mation,  which  usually  arise  from  intelligence  sources,  visual 
sightings,  acoustic  sensors,  etc.,  in  conjuction  with  (1). 

One  candidate  performance  measure  is  J,  the  negative  log 
posterior  distribution  of  data  partionings  evaluated  at  that  par¬ 
titioning  characterizing  the  correlator-tracker ,  at  hand.  A  closely 
related  measure  is  3',  the  negative  log  conditional  distribu¬ 
tion  of  the  observed  data,  conditioned  on  the  partitioning.  J' 
is  shown  to  have  a  much  simpler  form  than  J.  If  no  prior  know¬ 
ledge  is  available  conerning  what  the  true  partioning  of  data  is, 
then  a  uniform  distribution  assumption  lends  to  J  and  J-  dif- 
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fering  by  a  function  of  the  data  only  -  and  not  of  any  partitionings. 
In  this  case,  all  of  the  theoretical  justifications  for  J  equally 
apply  to  J' .  The  bulk  of  this  report  concerns  the  evaluation 
and  presentation  of  the  properties  of  3' . 

3'  can  be  expressed  as  the  sum  of  three  submeasures; 

Evaluation  of  J'"  involves  computing  Kalman  filters  for  the 
various  track  sets  of  the  correlator- tracker,  resulting  in  a  sub¬ 
measure  of  goodness  of  fit  of  the  tracks.  Another  submeasure 
.  computes  the  plausibility  of  the  false  alarm  set  determined  by  the 
correlator-tracker.  A  third  submeasure  determines  the  goodness  of  fit 

of  the  target  attribute  data  at  hand  (such  as  observed  ship  identi¬ 
fications,  frequencies  received,  hull  lengths  observed,  etc.). 

These  submeasures,  from  their  relatively  simple  structures,  are 
expected  to  be  implementable  for  real  time  incoming  data,  recur¬ 
sively  in  many  instances. 

In  order  to  properly  determine  3',  Kalman  filter  outputs 
of  the  correlator- tracker  being  evaluated  may  have  to  be  disre¬ 
garded  and  replaced  by  the  specially  structured  Kalman  filters 
used  for  3'.  This,  in  general,  is  due  to  either  or  both  of  the  fol¬ 
lowing  two  factors  being  present: 

1)  The  correlator- tracker  at  hand  may  in  its  own  opera¬ 
tional  form  be  using  an  adaptive  extended  -  or  in  some  way  non 
linear-Gauss-Markov  target  motion  and  measurement  models.  In  this 
case,  present  or  relatively  recent  target  positions,  e.g.,  will 

be  more  accurately  estimated  than  previous  ones,  and  correlation 
errors  can  actually  be  worse  than  for  non-adaptive  filters.  In  any 
case,  sums  of  the  apparent  quadratic  forms  of  data  fit  do  not  really 
represent  the  fit  of  the  entire  track  set  to  the  data. 

2)  Targets  may  be  maneuvering  so  that  the  linear-Gauss-Markov 
target  motion  assumptions  become  invalid. 

Appendix  C  exhibits  detailed  flow  charts  for  computations  re- 
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BACKGROUND 

In  brief r  the  correlation  problem  considered  here  may  be  de¬ 
fined  as  follows:  Determine!  within  feasible  computational  bounds, 
thab  partitioning  of  the  observed  data  set,  which  most  accurately 
reflects  through  its  component  track  sets,  the  true  (but  unknown 
to  the  observer)  partitioning.  In  the  correlation  problem  consi¬ 
dered  here,  data  is  accumulated  in  the  form  of  vector  valued  reports, 
labeled  as  to  sensor  system  source  and  time,  accumulated  over  some 
specified  sampling  interval.  The  component  sets  of  the  partition¬ 
ing  represent  either  track  sets,  each  corresponding  to  a  single 
target  source  (estimated,  if  the  partitioning  itself  is  an  estimate 
of  the  true  partitioning) ,  or  the  false  alarm  (or  clutter)  set. 

(The  latter  may  be  vacuous.)  The  observed  data  vectors  may  be 
classified  into  two  categories:  geolocational  and  non-geolocational 
attributes.  For  any  given  time  and  sensor  system  it  is  possible 
that  either  of  these  types  of  data  vectors  may  be  vacuous.  For  the 
geolocational  case  this  may  be  due  to  missed  detections;  however, 
generally  when  geolocation  information  is  not  missing,  the  number 
of  components  in  the  data  vector  will  be  fixed  and  typically  in¬ 
cludes  positional  and/or  velocity  observations.  False  alarms  may 
be  mixed  into  this  data.  For  the  non-geolocation  case,  the  number 
of  components  and  type  of  information  available  in  a  given  data 
vector  is  random.  Typically,  non-geolocation  target  attributes  can 
include  any  subset  of  the  following:  sighting  of  ship  identifica¬ 
tion,  observed  ship  silhouette/shape,  estimated  hull  length,  ob¬ 
served  flag  color,  number  of  sensors  of  a  given  type  detected  to 
be  on-board  vessel,  and  ship  machinery  frequency  outputs  observed. 

Broadly  speaking,  correlator- trackers,  i.e.,  schemes  which 
first  establish  partitions  of  the  data  and  then  use  these,  in  turn. 
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to  establish  tracks  of  the  targets#  nxay  be  divided  into  two  classes: 

•soft  decision*  and  'hard  decision*  types.  The  former#  basically 
establish  posterior  distributions  or  sets  of  possible  partitions 
to  be  considered#  generate  weights  based  on  the  likelihoods  or 
similar  measures  of  the  possible  partitions#  and  then  estimate  target 
motion  parameters  through  weighted  sums  of  Kalman  filters.  (See# 
e.g.#  Ref.  [2]).  On  the  other  hand#  hard  decision  correlator- 
trackers  may  be  identified  with  a  single  partitioning  of  the  data; 

■  the  estimated  target  motion  parameters  being  determined  by  a  Kalman 

filter  for  each  component  track  set  of  the  partitioning.  (See  Ref.  [3].) 

In  this  report#  only  those  correlator- tracker  schemes  are 

considered  which  can  be  identified  with  a  single  partitioning  of 
the  accumulated  data#  for  each  sample  time.  Thus#  this  includes 
all  hard  decision  correlator-trackers  and  those  soft  decision  ones 
which  can  be  suitably  modified,  such  as  by  choosing  that  partition¬ 
ing  of  the  data  which  maximizes  the  posterior  distribution  of 
possible  partitionings,  or  some  related  statistic. 

Extensive  overviews  and  general  analyses  of  the  multiple  tar¬ 
get  multiple  sensor  correlation  problem  may  be  found  in  NRL's 
Correlation  Handbook  [4] ;  see  also  the  earlier  work  of  Kullback  and 
Owens  [5] .  General  models  for  the  problem  are  presented  in  Sittler 
[6] ,  Reid  [2]  and  Goodman  [7] .  Bar  Shalom  [8]  and  Goodman  [7] 

Chapter  12  also  contain  technically  detailed  surveys  of  various 
correlator-tracker  schemes,  with  a  number  of  examples  presented 
for  both  the  hard  and  soft  decision  cases. 

STATEMENT  OF  THE  PROBLEM  -  PERFORMANCE  MEASURES 

Consider  the  problem  of  determining  how  well  a  particular 
correlator- tracker  technique  works.  It  is  obvious  that  if  target 
density  or  false  alarm  rates  are  sufficiently  high,  even  with  very 
accurate  sensing  systems  present#  no  correlator- tracker  will 
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perform  well.  Other  factors  also  contribute  to  correlation  diffi¬ 
culties  ^  such  as  sensor  reliability,  data  rates  and  variable  types 
of  target  motion,  including  xnaneuvers.  Despite  the  above  negative 
remarks,  there  appear  to  be  many  occasions  (moderate  data  rate 
with  few  targets  and  false  alarms  present,  e.g.)  when  the  correlation 
problem  may  be  amenable  to  reasonably  accurate  solutions.  In 
these  situations,  an  analytically  based  performance  measure  can 
serve  as  a  valuable  tool  for  ranking  and  comparing  quantitatively 
different  approaches. 

Due  to  the  possible  combinations  of  many  factors  -  including 
target  initialization  and  termination  times  within  the  sensor  areas 
of  interest,  detection  failures,  variable  target  motions  and 
maneuverings  and  matching  of  non-geolocation  attribute  information 
to  the  correct  track  sets  -  a  large  number  of  possible  partition¬ 
ings  of  the  data  into  track  sets  can  occur.  This,  in  turn,  can 
lead  to  an  exponential  growth  -  with  respect  to  sampling  time  -  of 
the  number  of  computations  required  in  implementing  or  monitoring 
error  bound  performance  of  correlator-trackers,  even  under  very 
simplistic  model  assumptions,  (See,  e.g.  Theorem  4  of  the  abridged 
version  of  Ref.  [7]  for  an  example  of  exponential  growth  in  the 
number  computations  needed  to  compute  target  motion  parameters. 

See  also  Ref.  [1].)  Thus,  alternative  performance  measures  must 
be  sought  which  do  not  involve  such  a  large  amount  of  computations. 

Performance  measures,  or,  as  they  are  often  called  MOE's 
(measures  of  effectiveness),  can  be  roughly  divided  into  two  basic 
classes:  (1)  observer  dependent  only  and  C2)  simulator-observer 

dependent. 
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A  performance  measure  is  defined  to  be  observer-dependent  only, 
if  it  is  a  true  statistic,  i.e.,  known  function,  of  the  cumulative 
data  and  the  partitioning  determined  by  the  correlator- tracker  in 
question. 

On  the  other  hand,  a  performance  measure  may  depend  in  imple¬ 
mentation  on  knowledge  of  what  data  constitutes  the  true  track  sets 
and/or  false  alarms,  i.e.,  on  the  true  partitioning  of  the  data. 

The  measure  may  further  depend  on  knowledge  of  the  correct  target 
motion  parameter  vectors  such  as  positions,  and  correct  time  of 
initialization  and  temination  of  the  targets  from  the  region  being 
surveilled.  In  this  case,  the  performance  measure  is  said  to  be 
simulator-observer  dependent. 

In  general,  an  observer  dependent  only  performance  measure, 
or  briefly,  score,  is  evaluated  relative  to  its  statistical  dis¬ 
tribution  as  a  random  quantity  (based  on  the  random  data) ,  or  at 
least  relative  to  its  statistical  moments,  conditioned  on  ^  assumed 
true  situation  (i.e.,  true  partitioning  of  data  or  even  true  tar¬ 
get  motion  parameters).  These  statistical  moments,  are  thus 
simulator-observer  performance  measures  generated  by  the  score. 

Some  examples  of  scores  are:  average  observed  track  depth 
or  persistence;  nximber  of  tracks  exceeding  a  depth  of  three  sample 
times;  and  niiraber  of  tentative  tracks,  i.e.,  track  sets  having  no 
more  than  a  prescribed  upper  bound  on  track  depth  (e.g.,  two). 

(See,  e.g..  References  [9]  and  [10].)  Another  measure  that  is 
observer  dependent  only  is  consistency  of  performance  (see  [10]). 

The  latter  measure  describes  how  much  difference  exists  between 
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the  initial  filtered-partitioned  output  data  from  a  given 
correlator- tracker  and  a  second  output. resulting  from  recycling 
the  Initial  output  through  the  same  correlator- tracker  (by  using 
it  formally  as  raw  input  data).  (See  also  Ref.  [1]  pages  6  and 
7  for  further  mathematical  details.) 

Simulator-observer  dependent  measures  of  performance  are 
much  more  plentiful  in  the  literature  of  correlator-trackers  and 
related  algorithms  than  those  in  the  former  class.  Wiener  et  al. 

[11]  present  a  number  of  examples  of  these  measures,  as  well  as  a 
general  overview  of  the  rule  MOE's  in  general  play  in  the  command 
and  control  aspects  of  ocean  surveillance.  (See  also  the  overviews 
in  [12]  and  [13].)  Typical  measures  here  include  all  measures  of 
accuracy  of  target  motion  parameter  estimators,  thus,  e.g.,  esti¬ 
mated  coverage  probabilities  of  target  parameters;  track  purity, 
i.e.,  percentage  of  correct  elements  in  any  given  track  set;  prob¬ 
ability  of  correct  association  between  new  reports  and  previous 
established  tracks,  and  percentage  of  targets  having  track  depths 
of  given  lengths.  Most  of  the  above  measures  may  be  averaged  with 
respect  to  sampling  times  and  nvimber  of  targets  present. 

This  section  is  concluded  with  a  brief  overview  of  the  litera¬ 
ture  of  simulator-observer  dependent  MOE's. 

Among  the  measures  Willman  [9]  uses  for  evaluating  the  per¬ 
formance  of  his  hard  decision  'matrix  scan'  correlator- tracker  are 
four  involving  actual  and  estimated  track  depths  (the  latter  obtained 
through  estimated  probabilities  of  consecutive  track  linkages) . 

In  addition,  Willman  employs:  a  measure  which  is  formally  the  same 
as  the  determinant  of  a  corresponding  two-by-two  table  of  entries; 
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percentage  of  old  targets  having  established  tracks  lacking  new 
assignment  of  reports;  and  percentage  of  misassigned  reports. 

Kullback  and  Owens  {[5],  pp.  31-42)  consider  numerical  dif¬ 
ferences  between  the  number  of  targets  and  false  alarms  estimated 
and  actually  present,  as  well  as  track  purity.  In  addition,  use 
of  non-geolocation  attributes  by  a  correlator-tracker  is  proposed 
to  be  measured  by  an  ambiguity  function  which  counts  the  total 
number  of  possible  ships  in  the  area  which  could  fit  those  attri¬ 
butes  assigned  to  each  track  set.  Other  measures  are  discussed  which 
measure  the  infommation  level  of  correlator-trackers. 

Adler  et  al.  [14]  have  proposed  a  number  of  interesting  MOE's, 
including  the  rate  of  track  fragments  (interruptions  in  tracking 
the  same  targets)  and  average  time  between  correct  correlations  (or 
linkages)  of  old  tracks  to  newly  arriving  reports.  Turner  and 
Harder  [15],  in  addition  to  the  usual  measures,  use  as  a  measure 
the  number  of  observations  required  of  a  target  to  form  a  firm 
track.  Reid  [2]  proposes  (in  addition  to  percentage  of  targets 
tracked,  lost,  etc.)  the  average  number  of  partitionings  of  data 
kept,  following  pruning  and  merging,  according  to  (his  soft  decision) 
correlator- tracker .  (See  also  the  MOE's  presented  in  References 
[16]  -  [18].)  The  performances  of  the  large  scale  correlator- 
tracker  presented  in  [19]  (resulting  in  various  model  improvements  - 
Personal  Correspondence)  were  measured  by  a  number  of  MOE's  pre¬ 
viously  mentioned,  including  track  purity,  accuracies  of  estimators, 
and  the  number  of  tracks  sets  having  relatively  short  depths. 

Observer  dependent  measures  are  presented  in  greater  detail 
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in  the  next  section. 

OBSERVER  DEPENDENT  PERFORMANCE  MEASURES  -  SCORES 

As  mentioned  in  the  previous  section,  many  measures  of  per¬ 
formance  exist,  but  in  effect  they  form  a  patchwork  quilt  of 
descriptions  for  the  behavior  of  correlator-trackers  in  general. 

Alspach  and  Lobbia  [20]  construct  an  observer  dependent  only 
performance  measure,  which  operates  on  a  correlator-tracker *s 
perceived  data  partitioning  only  through  the  total  number  of  re¬ 
ports  assigned  to  track  sets  and  the  total  number  of  reports 
assigned  to  the  false  alarm  set.  The  statistical  expectation  of 
the  evaluation  of  this  measure  (with  respect  to  the  randomness  of 
the  data)  holding  the  number  of  track  reports  and  false  alarms 
fixed,  is  also  a  performance  measure  which  is  both  observer  and 
simulator  dependent.  This  is  shown  to  be  a  concave  unimodal 
function  (under  certain  simplifying  assumptions)  of  the  total  nvira- 
ber  of  reports  assigned  to  track  sets.  This  function  is  seen  to 
be  a  linear  function  when  the  number  of  reports  assigned  to  track 
sets  is  less  than  or  equal  to  the  true  number  of  reports  that 
should  be  assigned  (via  the  true  partitioning)  to  the  track  sets, 
and  possesses  an  absolute  minimum  when  the  number  of  reports  decided 
to  be  targets  is  equal  to  the  number  of  actual  targets. 

In  this  report  a  score  J'  is  proposed  which  is  expected  to 
be  relatively  feasible  to  implement:  (twice)  the  negative  logarithm  of 
the  conditional  probability  function  of  the  possible  data  parti¬ 
tionings,  evaluated  at  that  partitioning  determining  the  correlator- 
tracker  to  be  evaluated.  An  associated  simulator-observer  dependent 
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performance  measure  (I')  is  also  presented.  This  measure  is 
essentially  the  statistical  expectation  of  J'  conditioned  on  any 
given  true  partitioning  of  data,  a  cross  entropy  measure. 

Although  Alspach  and  Lobbia  [20]  briefly  discuss  the  possible 
use  of  a  statistic  related  to  the  one  proposed  here  (but  do  not 
consider  non-geolocation  attributes) ,  they  dismiss  use  of  it  be¬ 
cause  of  apparent  difficulties  in  determining  the  non-random  terms 
consisting  of  determinants  of  innovation  covariance  matrices.  In 
effect,  in  Alspach  and  Lobbia 's  score,  the  non-random  terms  are  re¬ 
placed  by  the  product  of  the  number  of  data  points  decided  to  be 
clutter  or  false  alarms  employing  undetermined  constant  scores  for 
any  single  point  decided  to  be  clutter. 

In  Reference  [21]  and  [3]  (the  former  treating  data  recursively 
in  time,  the  latter  handling  data  non-recursively) ,  new  reports  are 
assigned  to  those  track  sets  and  to  the  false  alarm  set  such  that 
performance  measures  which  are  modifications  of  the  log  likelihood 
of  possible  partitions  -  a  statistic  related  in  form  to  that  pro¬ 
posed  in  this  report  -  are  maximized  (or  minimized) .  (See  Ref, 

[8J  for  concise  descriptions  of  other  correlation  techniques  which, 
in  effect,  also  use  performance  measures  for  determining  imple¬ 
mentations  . ) 

The  motivation  of  the  choice  of  the  negative  log  conditional 
or  posterior  probability  functions  of  possible  partitionings 
as  a  measure  of  a  correlator-tracker's  performance  is  based  on 
five  desirable  properties: 

1.  The  measure  in  its  initial  form  (before  taking  expectations 
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with  respect  to  the  data)  is  truly  observer-dependent  only,  and 
may  thus  be  used  as  a  real  world  scoring  method. 


2.  The  score  is  relatively  simple  in  form,  or  can  be 
reasonably  approximated  by  a  simple  structure,  suitable  for  real¬ 
time  implementations. 

3.  The  score  directly  reflects  the  goodness-rof-f it  of  the 
partitioning  of  the  correlator-tracker  in  question  to  the  given 
data,  employing  geolocation  terms,  non-geolocation  attributes  and 
false  alarm  data. 

4.  The  statistical  distribution  of  the  score  is  related  to 
a  chi-square  random  variable. 

5.  Use  of  the  score  can  be  justified  from  information  theory 
and  statistical  decision  theory  viewpoints. 

The  statistical  expectation  of  the  score  is  observer- simulator 
dependent  and  can  also  be  used  as  a  measure  of  performance.  In 
particular,  the  expectation  of  the  measure  will  have  smaller  values, 
generally,  for  correlator- trackers  which  use  more  information  and/or 
have  lower  risks  (equivalently,  better  approximate  the  Bayes  or 
maximal  posterior  data  partitioning) . 

The  above  five  properties  may  well  serve  as  a  general  guide¬ 
line  for  establishing  performance  measures  for  correlator- trackers. 
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OUTLINE  OF  THE  REPORT  . 

The  introductory  sections  describe  the  correlation  problem, 
in  general,  and  performance  measures,  in  particular.  'Hard'  vs. 

'soft'  types  of  correlator- trackers  are  detailed;  the  former  being 
more  conducive  for  being  evaluated  by  the  performance  measures  J 
and  J"  presented  in  this  report.  J  is  essentially  the  negative 
log  posterior  distribution  function  evaluated  at  that  partitioning 
Q  of  data  characterizing  the  given  correlator-tracker,  while  J' 
is  a  related  more  computationally  simpler  measure  involving  the 
conditional  distribution  of  the  data  (conditioned  on  Q) .  Surveys 
of  the  available  literature  for  both  observer  dependent  (J  and  J' 
are  in  this  class)  and  simulator-dependent  performance  measures  are 
presented.  Relations  between  these  measures  and  the  ones  proposed 
in  this  report,  where  applicable,  are  described.  A  guideline  is 
presented  for  desirable  properties  that  any  performance  measure 
for  correlator-trackers  should  possess.  (J'  satisfies  these  criteria.) 

The  Analysis  Section  first  establishes  (subsection  2)  com¬ 
prehensive  definitions  and  mathematical  models  for  all  assumed  tar¬ 
get  motions  (eq.  ( 2 . 1 )),  occurrence  and  location  of  false  alarms 
(subsections  2 ((3)),  ((6))),  and  observation  measurements  (eq.  (2.2)),. 
Following  this,  statistical  decision  theory  and  information  theory 
bases  are  established  for  use  of  the  proposed  scores  (subsection  3) . 

In  essence,  it  is  shown  that  the  lower  the  value  of  J  or  J',  the 
closer  the  correlator-tracker  (through  its  partitioning  of  the  data) 
matches  the  true  but  unknown  partitioning  of  the  data  into  the  cor¬ 
rect  track  sets  and  false  alarm  set. 

The  next  subsection  (4)  of  the  Analysis  Section  develops  the 
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full  structure  of  the  simpler  performance  measure  J'.  This  leads 
to  a  decomposition  of  J'  into  a  sum  of  three  terms,  reflecting 
goodness-of-data-fit  to  the  partitioning  determined  by  a  given 
correlator-tracker  (eq.  (4.1)).  The  geolocation  target  data  term 
is  shovm  to  be  the  sum  of  constant  and  quadratic  forms  of  innova¬ 
tions  that  are  outputs  from  Kalman  filters  (eqs.  (4.2)  -  (4.51)1. 

The  (geolocation)  false  alarm  term  is  the  sum  of  constant  and 
quadratic  forms  of  data  (eq.  (4.52)  -  (4.61)).  Both  terms  are  seen 
to  depend  explicitly  on  the  niimber  of  targets  and  false  alarms 
perceived  by  the  correlator- tracker ,  The  non-geolocation  target 
attribute  term  involves  discrete  sums  of  probabilities  (eqs.  (4.62)  - 
(4.70)).  If  a  normal  approximation  is  made  in  the  modeling  of 
these  attributes  (see  remarks  following  eq.  (4.70))  then  a  Kalman 
filter  (eqs.  (4.114)  -  (4.127))  may  be  used  for  evaluation  of  this 
term:  a  constant  plus  quadratic  forms  of  the  data.  (See  eqs. 

(4.63),  (4.64),  (4.68),  (4.69),  (4.79)  -  (4.82) .)  A  further  simpli¬ 
fying  approximation  -  which  avoids  basically  prior  knowledge  of  the 
randomness  of  the  non-geolocation  characteristics  is  given  in 
Appendix  B;  in  addition  the  accuracy  of  this  approximation  is  also 
demonstrated . 

The  final  subsection  (5)  of  the  Analysis  Section  develops  the 
distributional  properties  of  J',  It  is  shown  that  J",  conditioned 
on  the  partitioning  determining  the  correlator-tracker  being  evaluated 
is  distributed  as  the  statistical  independent  sum  of  a  constant, 
a  chi-squared  random  variable,  and  a  discrete  random  variable  (see 
eqs.  (5.1)  -  (5.5)).  Under  normal  distributional  approximations 
for  the  non-geolocation  attributes,  the  above  discrete  random 


variable  is  replaced  by  the  sum  of  a  constant  and  a  chi- squared 
random  variable  (see  remarks  following  eq.  (5.5).)  It  follows 
in  this  case  that  J'  is  distributed  as  the  sum  ol  a  constant  and 
a  chi-squared  random  variable  (eqs.  (5.6)  -  (5.9)).  It  is  also 
shown  that  J'  conditioned  on  a  partitioning  of  data  not  coinciding 
with  the  one  determining  the  correlator-tracker  being  evaluated, 
has  a  distribution  which  is  not  easily  computable  (not  even  a  non¬ 
central  chi-square  distribution) .  (See  equations  and  remarks  fol¬ 
lowing  eq.  (5.9).) 

Appendix  A  presents  a  procedure  for  calculating  matrix  in¬ 
verses  in  prescribed  block  form.  This  can  be  useful  in  evaluating 
parts  of  the  Kalman  filters  used  in  the  computations  for  the  geo¬ 
location  targets  and  non-geolocation  attributes  terms  (assuming  a 
normal  approximation) . 

A  further  approximation  to  a  quadratic  form  arising  from 
normal  distributional  approximations  for  the  non-geolocation  attri¬ 
butes  is  presented  in  Appendix  B,  where  also  error  bounds  are 
derived.  Use  of  this  approximation  minimizes  required  prior  know¬ 
ledge  of  the  randomness  of  the  true  non-geolocation  attributes  of 
the  targets. 

A  complete  set  of  flow  charts  for  computing  J'  relative  to 
incoming  data  is  presented  in  Appendix  C. 
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1.  INTRODUCTION 

In  order  to  be  able  to  define  and  analyze  consistantly  the 
performance  measures  proposed  here,  a  rigorous  mathematical  **  logical 
model  is  established  for  the  general  correlation- tracking  problem. 

(This  largely  simplifies  the  model  proposed  in  Ref.  171.)  This 
model  consists  of  eight  key  aspects:  ((1)}  sensor  systems,  ((2)) 

target  initializations  and  terminations,  ((3))  existence  of  target 
state  parameter  vectors  and  detection,  ((4))  partitionings  of 
observed  data,  ((5))  target  state  parameter  vectors,  ((6)}  observed 
geolocation  data,  ((7))  non- geo  location  attribute  data,  and  ((8)) 
total  observed  data. 

Although  the  model  obviously  simplifies  the  real-world  situa¬ 
tion,  it  is  expected  that  its  constituent  assumptions  are 
'reasonable*  approximations  of  reality.  In  keeping  with  the  attempt 
to  be  as  faithful  as  possible  to  the  real  world  setting,  the  model 
established  represents  a  relatively  short  data  sampling  period, 
j^,  and  can  be  changed  with  respect  to  each  new  sampling  period; 
all  of  the  available  output  information  from  the  previous  sampling 
period,  being  used  as  input  -  prior  information  for  the  next  period. 
This  disjointing  of  the  sampling  times  into  short  segments  should 
make  more  valid  the  homogeneous  linear  Gauss-Markov  motion  and  measure¬ 
ment  models  used  in  ((5))  and  ((6))  .  The  latter  assumption  can 
account  to  some  degree  for  variations  in  types  of  target  motion  by 
allowing  in  the  model  a  reasonably  large  state  vector  dimension 
with  zeros  possibly  occurring  in  particular  entries.  (Thus, 
quadratic  polynomial  motion  includes  as  a  special  case,  straight 
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line  motion,  with  the  possibility  of  estimated  error  covariances 
being  actually  larger  than  necessary.)  This  results  in  the  avoid¬ 
ance  of  a  potentially  large  branching  problem  that  can  arise  in 
attempting  to  model  the  general  correlation-tracking  problem. 

(See,  e.g.  Goodman  [7J ,  especially  Theorems  4  and  5  of  the  abridged 
version  for  an  illustration  of  the  complexity  arising  when  variable 
target  motions  and  maneuvers  are  modeled.) 

Even  over  a  relatively  short  -sampling  period,  a  target  may 
engage  in  maneuvering  -  such  as  zig-zags  or  circular  motion  -  which 
is  not  really  modeled  by  the  same  linear  Gauss-Markov  target  motion 
model.  Yet,  a  given  correlator-tracker  may  still  retain  the  ability 
to  follow  that  target,  and  thus  essentially  put  all  observations 
of  it  into  the  same  track  set.  (This  will  usually  be  carried  out 
by  highly  nonlinear  adaptive  Kalman  filters  -  which  do  not  reflect 
the  target  motion  model  assumed  here.)  Clearly,  a  linear  Gauss- 
Markov  fit  to  this  track  set  -  which  is  one  of  the  computations 
required  to  obtain  the  score  of  the  correlator-tracker  in  question 
(see  eq.  (4.20)),  -  in  general  for  this  situation  is  not 
appropriate . 

Consequently,  a  special  procedure  is  used  in  computing 
the  score  for  a  correlator-tracker  which  has  at  least  some  track 
sets  generated  apparently  by  a  maneuvering  target  as  described 
above.  In  essence,  the  procedure  retains  the  mathematical  rigor 
of  the  model,  by  simply  replacing  the  single  track  set  in  question 
by  a  disjoint  union  of  different  (but  almost  contiguous)  track 
subsets,  each  distinct  track  subset  based  on  measurements  eminating 
from  a  different  linear  Gauss-Markov  motion  model  (as  is  permitted 
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in  the  overall  model) .  The  termination  time  of  each  such  subset 
(and  consequently,  the  next  sampling  time  being  the  initial  tiige 
of  the  following  subs’et)  is  determined  by  either  a  simple  chi- 
squared  test  involving  goodness  of  fits  or,  more  generally  by  using 
the  non-geolocation  attribute  of  common  -  but  unknovm  -  identity 
in  monitoring  the  entire  goodness  of  fit  data  probability  function. 

See  Appendix  C,  for  implementation  of  this  procedure  and  for 
overall  flow  charts  for  computing*  the  score  J'  of  a  given 
correlator-tracker  as  a  function  of  incoming  data. 

It  should  also  be  noted  that  often  correlator-trackers  operate 
with  nonlinear  adaptive-extended  Kalman  filters  of  the  track  sets. 
In  these  cases  the  original  filters  must  be  replaced  by  the  linear 
ones  for  the  model  developed  here  (see  eqs.  (2.1),  (2,2))  in  order 
to  reflect  the  total  goodness  of  fit  of  the  track  sets. 
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2.  ASSUMPTIONS  AND  BASIC  DEFINITIONS  FOR  THE  CORRELATION-TRACKING  MODEL 

(( 1 ))  Sensor  Systems 

q  known  fixed  sensor  systems  sample  data  in  the  form  of  vec¬ 
tor  reports  at  known  and  possibly  random  outcome  times 
t^ <  t,  <  . . .  <  t .  ,  where  j.  is  chosen  relatively  small  based  on 

u  3o  0 

experience  or  other  factors.  (See  the  discussion  on  choice  of 
in  the. previous  :  subsection.) 

((2))  Target  Initializations  and  Terminations  in  API 

An  unknown  ntunber  M^^^  (to  the  observer)  of  targets  begin 
existence  in  the  sensor  areas  of  interest  (AOI)  to  time  t^. 

This  is  due  to  targets  either  entering  the  AOI  for  the  first  time 
and/or  becoming  sufficiently  active  that  they  may  be  detected  by 
at  least  one  of  the  sensor  systems  present.  Some  of  the  targets 
in  the  AOI,  later  may  terminate  existence  (relative  to  the  AOI) , 
i.e.  exit  the  AOI  or  quiet  down  so  that  they  can  no  longer  be  de¬ 
tected  by  any  of  the  q  sensor  systems  present.  The  number  of  tar¬ 
gets  actually  existing  ^  t^  is  denoted  by  M^,  also  unknown. 

Thus,  Mj^M^^\  M^^^  is  non-decreasing  in  j,  and  and  M^^^ 

are  random  integer  outcomes. 

Related  to  the  above  definitions,  associated  with  any  target 


i  s  M<'>  il  0  M<J> 

j=l 


are  two  unknown  random  integers  u^,  v^. 


0<u.<v.,  where  t  is  the  time  of  target  i's  initial  existence 

'^i 

in  the  AOI  and  t^  is  the  time  of  its  termination.  It  is  assumed 

that  any  target  exists  over  successive  sampling  times  and  if  it  ceases 
existence  (for  two  or  more  sampling  times)  and  begins  again  at  a 
later  time,  it  is  considered  a  distinct  target  here. 
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i(3)J  Existence  of  Target  State  Parameter  Vectors  and  of 
Observed  Data  Vectors  (Detection) 

For  any  target  i,  time  t^  and  sensor  k,  represents 

its  state  parameter  vector  (not  dependent  on  k) , 

observed  geolocation  measurement  vector,  and  Y-  its  observed 

non-geolocation  attribute  measurement  vector.  ^Qjj^  set  of 

false  alarm  vectors  occurring  at  time  t^  due  to  sensor  k. 

X .  .  iff  target  i  exists  at  time  t . ,  in  which  case 

dim(X..)=  m  is  known.  Thus  X..J^^  iff  u.<  j<v. 

The  following  is  closely  connected  with  sensor  detection: 

Z...  ^  ^  iff  sensor  k  makes  at  time  t.  one  geolocation 

ijK  ^  3 

measurement  of  target  i,  in  which  case  dim(Z^jj^)  =  r^j^  is  known. 

^Ojk  sensor  k  receives  at  time  t ^ ,  f^j^  ^^jk  = 

geolocation  measurements  of  false  alarms  in  which  case,  the  false 


alarm  set  is  H  |  a  =  1,2, . . . ,  f  ,  where  Z^^^^  is 

the  false  alarm  vector  seen  by  sensor  k  at  time  t.  and 

dim(ZQjj^^)  =  r^j^  is  known,  f is  an  unknown  random  outcome, 
distributed  exponentially  as  Expo  (X..),  X.,  known,  f.^'s  are 

3X  3K  3K 

statistically  independent  with  respect  to  different  (j,k)*s. 

Y.  ..  ♦  iff  sensor  k  makes  at  time  t.  one  non-geolocation 

iDK  3 

target  attribute  measurement  of  target  i,  in  which  case 
1  ^  dim(Y^jj^)$  b,  dim(Y^jj^)  and  b  are  known. 


is  always  unobserved. 


K 


ID 


^ijk  always  observed. 

Index  i  ^0  (including  i  =  0  for  the  Zq^j^'s)  is  unknown 

to  the  observer,  while  j  and  k  are  always  known.  Index  u  is 
always  unknown.  However,  when  the  observer  forms  a  partitioning 
of  the  data  (see  the  next  subsection),  relative  to 
assxuaed  to  be  formally  true,  index  i  (including  i  =  0)  becomes 
(formally)  known. 

If  X.  .  =  (|»,  then  clearly  Z.  ,,  =-y.  =  (k ,  for  1<  k<  q. 

On  the  other  hand,  ^  (fr  does  not  guarantee  and/or 

being  nonvacuous.  For  example,  sensor  k  may  miss  a  dis¬ 
tinct  detection  due  to  resolution  problems,  noisy  background,  or 
reliability  problems  caused  by  equipment  or  human  operator  failure. 

(j).  . 


Define,  for  any  i;  l<i<M  -^  ;  j,  k, 
Z 


-  Wjk j  ' 


z  ^  z 

Ojk  —  ^Ojk 


Z..  represents  all  observed  data  of  target  i  (for  i^  1) 

1 JK 

at  time  t^  by  sensor  system  k.  It  is  called  a  report. 

((4))  Partitionings  of  Observed  Data 
For  any  i,  l^i^M^^\  define 


0^(3)  ^  { (i,  a ,k) I  for  all  a,  k, 

0^a4j  &  l<k$q  such  that 

^iak  ^ 

-  target  track  index  set  1  up  to 
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s=  (  (O#o,k,(i))  I  for  all  a,  k,  u. 


0^a$j  &  lgk<q  &  l^w  ^ 
Vk^*  (and  hence  f„j,>l)> 


&  such  that 


s  false  alarm  index,  set  up  to  t^  (possibly  vacuous) . 
The  partitioning  index  of  data  up  to  t^  is  given  by 


Q<j)  «  ,Q.«)|  ieA<l>)  , 


where 


Mz  ii 

corresponds  to  the  set  of  all  distinct  target  track  sets  established 
up  to  tj,  including  the  false  alarm  set  as  a  special  track  set  (0) 


Notice  that  Q 


is  determined  by  the  (detection)  set  of 


all  (i,a/k)'s  for  which  indirectly  by  all  T^j^^'s, 

by  all  u.'s,  v.'s;  and  by  the  set  of  all  (0,ark,(o)'s  for 

which  f  ,^1. 
ak 

Given  index  there  are  infinitely  many  corresponding 

partitioning  outcomes  Q^^^(Z^^^),  depending  on  the  values  of  the 

Z  .  .  '  s ,  where 
lak 

Q(j)  ,2<J),  df  |i,A<3>)  , 

Q^(J)(Z<J),  H  (z^  ^  I  (i,a,k) 


In  general,  we  identify,  if  ambiguity  does  not  arise, 
q(3)(Z<3>)  with 

Note  that  all  partitionings  are  unlabeled, that  is  the 

indices  i  for  each  track  set  (or  false  alarm  set,  for  i =  0) 

does  not  identify  what  target  i  really  is.  It  is  just  a 
convenient  index. 
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In  some  of  the  following  equations  equality  holds  only  after 
a  suitable  rearrangement  of  indices  is  made. 

((5))  Target  State  Parameter  Vectors 
The  results  here  do  not  depend  on  : 

Given  outcomes  (u.,  v.)  for  i  =  l,..,  and  equivalently 

those  (i,j)‘s  for  which  X..  ^ 

13 

~  I  u^  <  a^minCVj^,  j)  } 

==  ^^ia  l0^“^3  &  } 

=  set  of  state  vectors  of  target  i  up  to  t^. 

Xj  M  {X^j  I  l^ig^M^^^  &  x^^  5^  (J>} 

=  set  of  all  state  vectors  of  targets  existing  at  t . , 
X^^^  M  (X^  1  0^«  ^  j  &  X^  ifr> 

—  set  of  all  state  vectors  of  targets  existing  sometime 

up  to  t . . 

3 

Initial  state  vector  X.^^  of  target  i  (at  t^  )  is  assumed 

i  “i 

to  be  distributed  normally  N^(E(X^  Cov(X.,  )),  where 

Cov(X^  )  is  known  and  E(X.  )  is  unknown,  unless  otherwise 
specified. 

(Prior  information  for  3^u^,  may  be  used  to  determine 

both  or  either  of  its  moments , ) 

If  °  “ '^i+i' •  • '^i'  assume 


X, 


=  *X.  ^  +  G  'W. 

a  1,0-1  ^<x  ”ia 


(2.1) 


linear  Gauss-Markov  homogeneous  target  motion,  where 
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A  is  m  by  m  known  transition  matrix, 
o 

6^  is  m  by  n  known  coefficient  matrix. 

W.  is  distributed  N  (0,  P  );  P  known, 
xa  n  a  o 

W^^'s  are  statistically  independent  for  different  (i,a) 's  ; 
they  are  all  unobserved  random  variables, 
and  are  statistically  independent, 

for  different  i's  are  statistically  independent. 


{(6))  Observed  Geolocation'  Data 

Zi<3>  as  1  0<a  <j  s  l<k  <g  4  ♦),  for  any  i. 


iak 

0  <  i  <  , 


=  geolocation  data  for  track  set  i  (and  hence  for 


^ij 


Z . 
3 


target  i ,  i  unknown  to  the  observer)  up  to  time  t ^ , 


M  {Z.  .. 
ijk 


1<  k  <  q  4  Z^.^J‘  .}, 

th 


=  geolocation  data  for  the  i  track  set  at  t^ 
(from  all  sensors) , 


df 


(Z^jj^  I  0<  i  <  &  l^k  iq  & 


jk 


=  geolocation  data  for  all  track  sets  at  t ^ , 


2  (j)  ^  {z  I  0  <  i  ^ 

lak  ' 


4  OiaSj  4  l<k<q  4 


=  {Z^  1  O^a  <  j  &  Z^  5^  4.) 


=  geolocation  data  for  all  track  sets  up  to  t^ 

=  total  observed  geolocation  data  up  to  time  t ^ . 
Define  similarly,  etc. 

For  any  ^  the  observed  geolocation 

target  data  or  measurement  is  assumed  to  have  a  linear  regression 
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.  relationship  with  respect  to  a  corresponding  target  state  vector: 


"ijk  =  *  ''ijk 

In  more  compact  form, 


Z..=B.  X..+V..  ; 

13  i3  ID  ID 


df 


hi] 

Kn\ 

: 

,v. .  ~ 
ID 

for  l|Dc<q 
■  such  that 

• 

IV.  . 

’'ijk  ♦ 


for  l^k  < 
such  tha^ 


^ijk^  4‘. 


is  always  a  known  r^j^  by  m  matrix  *  '^ijk 
unobserved  and  distributed  statistically  independent  of 
indeed  of  X^^^;  statistically  independent  for 

different  (i,j,k)'s. 

^i-ik  distributed  normally  as  N  (0,  R..),  with 

J  Tjk  3k 

always  known. 

Thus,  for  i  =  l, 

V..  is  distributed  as  N  (0,  R.), 

^  J  S  •  •  •  T 

i3  1-^ 

where 


R.  M 


"ij 


such  that 
Zijk  ^  ^  r 

^  r • .  =  dim ( Z . . ) 


I  such  thatl 

\^ijk  ^  7 

Since  r^j^  is  known  to  the  observer »  if  is  given 


(2.2) 


and 


,  then 
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Note  that 

=  total  observed  geolocational  false  alarm  data  up  to 

^j* 

For  i=0  such  that  1=“  =  * 

it  is  assumed  that  is  statistically  independent  with  respect 

to  different  (j,k,u)'s  ,  all  other  l^i^M^^^  ,  and  to 

.  In  addition^  it  is  assumed  that  i  f^j^  $=1,  for  fixed  j,k, 
the  ^Qjj^^’s  for  1^  are  identically  distributed  as 

^Ojki  =  "‘Sjk'  V'  ®jk-®"^0jkl)  “jk  =  '"'"""ojkl' ' 

and  M.,  are  assiimed  known.  Although  in  general  f..  is 

3K  3K  3K 

unknown  to  the  observer,  given  (true)  f.,  becomes  known. 

3^ 

Similar  properties  are  assumed  for  the  ^ojko)'®  9ivei^ 

((7))  Non-qeolocat ion  Target  Attribute  Data 
■  Let  a  and  b  be  fixed  known  positive  integers  and 
C  =ss  {Cj^,  3  fixed  known  set  of  distinct  objects  where  each 

;  1<  a  <_  a. 


is.e  i  b 


The  C  a's  represent  possible  non-geolocation  attributes  of 

ap 


any  target. 
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For  example  for  b  =  3, 

C  ,  =  visual  identification 
al 

where  each  is  a  distinct  ship  name, 

C  ^  =  observed  flag  color 

e  {Red,  Green,  Brown,  White}  , 

of  radar  sensors  determined  to  be  on-board  the  ship 

£  {0,  1,  2,  ••,  10}  • 

Thus,  if  all  possible  combinations  of  attributes  can  occur, 
then  a  =  100* 4 *11  =  4400;  but  if  certain  combinations  of  values  of 
C  C  C  .  are  excluded, then  a  may  be  a  good  deal  less  than 

ct  X  Ok  z  a  3  ^ 

4400.  (This  will  often  be  the  case.) 

Define,  random  variable  (not  dependent  statistically  on 

over  the  set  C  to  represent  the  distribution  of  possible 
attributes  occurring  for  target  i,  i=l, 


The  prob¬ 


ability  function  of  is  known.  (H^'s  outcome  is  unknown  to 

the  observer.)  Also  assvune  that  the  distribution  function  for 
is  not  dependent  on  t ^ ,  for  the  duration  of  target  i ' s  existence 
in  the  A.O.I. 


Now  for  any  i,j,k  such  that  and  1^  k  ,  let 

T^jj^  be  always  an  observed  random  subset  of  {1,  ..,  b}  which  is 
statistically  independent  for  different  (i,j,k)'s  and  statistically 
independent  of  and  t«4- 


Let  a and 


arbitrary,  1  ^  o',  a"  ^a  and  let  Y'  —  (C  ,  „)„  _ 

a  ,6  e  e 

—  (C  .) 

“  '  ®  1<  e 

Define  conditional  random  variable  Y. by 

13K 
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prCY-.k  »  ®  =  P(Y"  ,  Y-*  j 


(2.3a) 


ijk 


a  knovm  function  of  its  arguments  (not  dependent  functionally  on 

i  nor  k  and  only  dependent  on  through  the  outcome  of 

T.  . 

X3k 

Analogous  to  the  definitions  in  subsections  ((5))  and  ((6)} 
define  Y..,  Y^^)^  etc. 

•V 

^ijk  considered  the  observation  of  H^,  the  actual 

non-geolocation  attribute  present  (given  selection  outcome  T. ..  ) 

13K 

of  target  i , 

Assume  that  given  the  all  of  the  (Y^^k  ) 's 

are  statistically  independent  with  respect  to  different  (j,k)*s; 
suppose  all  I  Q^^^)'s  are  statistically  independent  for 

different  i's» 


Equivalently,  for  outcome  T. 

IjK 


y. =  B. ..  »H.  +  V.  - 
13k  i3k  i  13k 


In  more  compact  foirm. 


Y. .  =  B. 'H.  +  V. . 
13  p  1  ID 


(2.3b) 


\  ijq/  for  <q  \  ijq  /  for  1  <  k  <  q 


such  thaF 

Y.  M-  * 
13k  ~ 


such  that 
^ijk  ^  ♦ 


R  df 

13k 


6  •  •  1 . 

i3kl 


8. 


fi  ■  -  (Of  • *  t  Of  If  Of  ..f  0)f 

“ijkB 
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is  always  known,  where  the  1  in  the  above  1  by  b  row  vector 
occurs  at  only  the  position;  for  outcome 


’'ijk  “  *’ijk;2'  ' 

“ijk  —  <’'ijk>' 

The  unobserved  discrete  random  variables  with 

known  distributions,  and  are  statistically  independent  of 
For  different  (ijk)  's,  the  are  statistically  independent 

Given  a  ir  Vnoi.,n . 


,  is  known. 


Unless  prior  information  is  available,  the  probability  dis¬ 
tribution  function  for  any  is  not  really  known  to  the  observer. 

Consequently  for  purposes  of  implementing  the  proposed  score  in  a 
real-world  setting  (where  the  target  identifications  are  not  known) , 
it  is  assumed  initially  that  all  of  the  H^*s  are  statistically 
independent  with  respect  to  different  i's  and  possess  identical 
uniform  prior  distributions  over  C.  On  the  other  hand,  for  a 
simulator-dependent  measure,  the  H^^'s  may  be  assigned  different 
known  prior  distributions  (including  possible  dirac  ones) ,  if  both 
the  attribute  properties  of  each  target  are  known  and  the  index  i 
is  identified  with  the  proper  target  by  the  simulator. 

Thus,  we  define  here 


pr(H.  =  y"*  ! 


=  pr(H^  =  Y  ) 


1/a,  if  Y''e  C 
0  ,  if  Y'jr  C 


(2.4) 
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'*\>  vjk  ■  V**** 

fits))  '  Totftl^'Obsaryod  Data  and  Additional  Notation 
By  aa8UBq|>tions  and  from  subsection  ((6))  it  follows  that 


prC^^^  I  .  (2.5) 

=  pr(Z^^^  I 
=  pr(Z^^^^  I 

noting  again  that  since  is  given,  so  are  the  outcomes  of  the 


X.  ,  ,  . 

xak's 

By  suitable  rearrangements,  the  total  observed  data  vector 
may  be  broken  up  into  geolocation  and  non-geolocation  atrribute  com¬ 
ponents,  and  further  into  geolocation  observed  target  data,  geo¬ 
location  observed  false  alarm  data  and  non-geolocation  attribute  data: 

.  (j)\ 

Z<:j^  =  =  i  (2.6) 


Using  the  Z  notation,  the  following  interpretations  hold: 

=  track  set  i  (i.e,,  all  data  corresponding  to 
up  to  tj, 

Z^jj^  =  report  (ijk) 

=  all  data  for  track  set  i  at  t^,  from  sensor  system  k, 

^ij  "  ^^ijk^l<  k  <  q, 

Z^  =  all  data  at  t^, 

Z^^^  =  all  data  up  to  t^,  etc. 
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Finally,  notation  will  be  introduced  which  emphasizes  at  each 
new  sampling  time  just  prior  to  forming  a  new  updated  partitioning 
(which  will  incorporate  and  be  based  upon  the  new  data) ,  the 
observer's  lack  of  knowledge t 


j  m 

k  =  1, , , /q 


C2.71 


up  to  suitable  rearrangments  of  2  where  m . .  i  0  is  the  total 

J 

number  of  data  reports  observed  by  sensor  k  at  t ^  -  always  known 

where,  recalling  from  subsection  ((3)) ,  a  single  report  consists  of 
either  (geolocationl  data  ^^jk  non-geolocation  attribute 


aata 


for  sensor  k  at  t^ . 


is  of  the  form 


H  t «  •  • « 

Y3k  13k 


(r^j^  by  1) 


(2.8a) 


for  some  corresponding  unknown  (to  the  observer)  i  =  i(Y,j,k), 


1  <  i  < 


(M 


(j) 


also  unknown) ,  or  it  is  of  the  form 


^Yjk  “  ^Ojka;  ^ojku)  ^^jk 


(2.8b) 


where  a  =  u(Y/jfk),  l^ui^f^j^,  the  total  number  of  false  alarms 
observed  by  sensor  system  k  at  t ^ .  (Recall  that 

=  ^^Ojko)  I  ^  =  “  =  ^jk^*  roost  one  report  B  corresponds 

to  each  target  for  sensor  k  at  t^  ;  the  remaining  reports  are 
false  alams.) 

Note  than  for  any  ,  m^j^  does  not  depend  on  , 

but  does  satisfy  the  relation 
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3.  BASIC  DECISION  AND  INFORMATION  THEORETIC  PROPERTIES  OF  THE 
PROPOSED  SCORE 

One  basic  score  J  for  a  given  cor relator- tracker  is  defined 
here  to  be 

j(Q(j)^  Z^j))  ^  -  2  log  pr(Q^3)  =  Qt3>|z^3)  =5<j))  (3.1) 


where 


6<J> 


is  that  outcome  of  partitioning  of  data  corresponding 


•  (-i) 

to  the  correlator-tracker  in  question  and  Z  ■*  i 
data  outcome. 


is  the  observed 


The  corresponding  simulator -observer  dependent  measure  is  de¬ 


fined  to  be 


(3.2) 


where  Q  (Z  •*  )  is  an  outcome  of  the  true  (but  unknown  to  the 


observer)  partitioning  of  the  data. 


Thus  I(  )  is  a  cross  entropy  measure.  (See^  e.g,.  Ref. 


[22] .) 

Consider  the  statistical  decision  theory  game  with  parameter 

0  (i) 

and  decision  space  being  the  set  of  all  possible  outcomes  Q 


of  Observed  data  Z^^^ 


has  distributions  of  all  relevant 


random  quantities  determined  from  the  previous  assiamptions .  For 
this  game,  the  loss  function  L  is  of  the  zero-one  type:  for 


If 

decision  outcome  and  true  parameter  value  Q 


of 


1  iff 
0  iff 


(3.3) 


Then  the  Bayes  decision  function  for  this  game  is  identifiable 
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iJJi  ^ .  - 


with  th«  posterior  naxinvmi  likelihood  estimator  ' 

it  •  O  •  p 


m.*  pr(Q<3>  .  Q<3>  I  Z'l'  =  5'’>) 

occurs  for 

(See  Ref.  [23],  Chapter  11  and  Ref.  [24]  for  background  and 
elaboration  of  results.) 

Thus  it  immediately  follows  that 


over  all 


5(j) 


(3.4) 


(j)  _ 


uniquely , 


min  J(Q 


) 


•  (i)  -  (i) 

occurs  for  Q'^-"  =  Q'-"  ,  uniquely. 

If  the  above  statistical  decision  game  is  modified 
between  two  hypotheses  (the  results  are  readily  extended 
than  two) 


(3.5) 

to  choosing 
to  more 


H 


(1) 


Q(j)  ^  5(j) 


vs.  H 


(2) 


Q(j)  =  'Q-(j) 


with  loss  function  and  prior  distribution  specified  by  (assuming 
Q<^>  Q<^>) 


L(Q<1>,  Q«i>)  « 

L(Q<1>,  Q«>)  2£  1,(11) 

L(Q<1>,  q(J>)  Ml'^D 
L(Q<J),  Q(i)) 

l(12),  l<21)  >  L<11),  l'“) 


(3.6) 
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(3.7) 


pr(Q<3)  =£(3), 


pr(Q 


(j)  =  ^ 


Q'-'M  -P2  =  1-p 


1  ' 


then  the  Bayes  decision  function  for  this  game  is 


Decide  H 


(1) 


pr(z'J>  ■.q"'!*) 


>  t; 


Decide  H 


(2) 


^  pr(z‘3>  =i‘3>|Q'3>  =5'‘5>) 


where  threshold 


df  P2  ..L<2P-Lt22) 

^0  —  ^  j.TrfT’7i7^T) 


(3.8) 


(3.9) 


Equivalently,  taking  logs  (see  the  definition  of  J'  in  eq. 
(3.14)) 


Decide  H 


(1) 


iff 


,Z  ^^^ )  -r  J' 2  log  T 


Decide  H 


(2) 


iff 


J.  (gCj)  Z^3)j  _  J- (q"^3)  JCji)  2  log  T'  . 


(.3.10) 


(For  pj^  =-p2  =  h  and  ~  ^21  ^11  ^22'  *^0  ~ 

^2  log  Tq  =  0.) 

Note  the  equivalent  decisions 
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Decide  iff 


<  J(Q*^^^  -  2  log  T 


Decide  H^2) 


-  2  log  T 


(3.11) 


where 


df  Lt21)_Lt22) 


L'  '  -  L' 


(3.12) 


For  a  discussion  of  the  distribution  of  the  J'''s  and  J's, 
and  hence  the  type  I  and  type  II  decision  errors ^  as  well  as  the 
probabilities  of  correct  decisions,  see  subsection  5. 

Suppose  now  the  prior  distribution  of  is  uniform  over 
the  set  of,  say,  Yj  possible  outcomes. 


Then 


cr(Q<3>,z<5>)  - 


where 


(3.13) 


«  .2109  pr<Z<l>  =  l'^>|0<3>  =  Q‘i>)  (3.14) 


Dj(z‘^’)'S  2  log  (Tf  -pr(Z^^’  =  Z*^’)) 


(3.15) 


5<J>. 


Thus  J  and  j"  differ  by  a  function  of  and  not 


Hence,  in  this  case,  defining  =  Q^^^(Z^^^)  as  the  condi¬ 


tional  maximxim  likelihood  estimator  of  ,  i.e., 


(over  all\ 

Q^>  ) 


maxpr(Z<3>  = 


(3.16) 


* (i)  ^ (i) 

occurs  for  Q  =  0  ,  uniquely, 
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then  equations  (3.5)  and  (3,13)  imply 


^  min  .  J(Q^^^  ,Z  and  min  J'(Q^^^,Z^^^) 

(over  all!  /over  all\ 

,  Qt:j>  /  A  6^3)/ 

occur  for  Q'J'  =  q'-*'  =  uniquely. 

In  addition,  eq.  (3.13)  implies  that 


where 


(3.17) 


(3.18) 


(3.19) 


2£E^^.jCD.(Z«Vlrf^>=6'^3V 


(3.20) 


the  latter  being  a  function  of  Q 
Then 


only. 


min  and  min  i ,Q ^^^ ) 

(over  all\  /over  all\ 

,  J'i’j  V  o'"’/ 


(3.21) 


" (i)  ' (i)  " (i) 

occur  uniquely  for  the  same  function  Q'-*'  =q'J'  =0  "^  ,  by  using 

equations  (3.17)  and  (3.18).  Also,  the  Fundeunental  Inequality  of 

Information  Theory  (see,  e.g..  Ref.  122],  Chapters  2  and  3)  implies 
•  (i ) 

that  if  the  Q  's  are  restricted  to  constants  (not  non-trivial 

functions  of  Z^^^),  in  the  two  equivalent  minimizations  in  equation 

(3.21),  the  corresponding  mimina  occur  for  the  common  value 
*  (i)  *‘’(i) 

Q\j/  _  Q'J',  These  minima  are  easily  seen  (via  eq.  (3.4) ,  e.g.) 
to  have  larger  values  in  general  than  those  for  the  unrestricted 
case  in  (3.21). 
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In  sxffiunary,  the  above  equations  imply  that  if  correlator- 

* (i)  * (i) 

tracker  has  partitioning  output  Q  operating  on  data  Z 

and  the  true  but  unknown  partitioning  of  the  data  is  actually  out- 

"Hi  * 

come  then  the  lower  the  score  J(Q  -'  ,  Z  -^  )  -  noting,  if 

the  uniform  prior  distribution  assumption  for  is  made,  J 

can  be  replaced  by  J'  -  the  better  Q'-*'  approximates  Q 
in  both  decision  and  information  theory  contexts. 
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4,  COMPUTABLE  STRUCTURE  OP  THE  SCORE 

Throughout  this  subsection  it  will  be  assumed  that  the  prior 
distribution  of  is  uniform  so  that;  without  loss  of  generality, 

we  consider  only  the  score  3'  and  its  expectation  (with  respect 
to  data  i".  For  ease  of  notation  here,  outcomes  of  random 

quantities  -  such  as  Q  with  respect  to  Q  -•  (partitioning) 

with  respect  to  Z^^^  (observed  data),  etc.  -  unless  ambiguities 
arise,  will  be  identified  with  the  corresponding  random  quantities. 

The  structure  of  3'  is  seen  to  be,  using  eq.  (2.5)  and 
the  calculus  of  conditional  and  joint  probabilities: 

j.  (q(  j)  ^2  ^^^ )  =  -^2  1og  pr  ) 

-  2  log  prCZp^^^-lQ^^^ 

-  2  log  pr(Y^3)  ,  (4.1) 


a  decomposition  into  a  sum  of  three  terms:  the  first  involving  geo¬ 
location  target  track  data,  the  second  pertaining  to  false  alarms, 
and  the  last  relating  to  the  non-geolocation  target  attribute  data. 

In  turn,  each  of  the  terms  in  eq.  (4.1)  may  be  decomposed  by 
straightforward  use  of  the  special  notation  developed  in  subsection 
2,  and  eqs.  (2.1)  -  (2.4),  with  the  corresponding  assiunptions : 
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is  the  measure  of  overall  goodness  of  fit  of  the  geolocation  data 
for  partitioning  to  all  the  track  sets  (determined  by 

(Zijj|Zi^“  represents  the  random  variable  of  goodness 

of  fit  of  the  geolocation  data  for  partitioning  to  track  set 

i  at  t^  (based  on  previous  data) ,  and 

^  -2  log  pr(Z^^lz^^“"^lQ^^^  (4.3) 

is  a  measure  of  goodness  of  fit  of  geolocation  data  at  t^^^  to  the 
i^^  track  set,  with  respect  to 

is  the  i^^  track  set  goodness  of  fit  of  the  geo¬ 
location  data  for  (corresponding  to  target  i,  the  index  i 

being  unknown  to  the  observer) ,  and  is  given  by 


For  a  =  0,  )  may  be  interpreted  as  either 

or  using  information  as  (Z^q  I  Z^  ,Q ^^^)^ 

If  Z .  ^  <|)  and  Z.  ^  4>  r  then  since  X.  ,  .V.  and 

lOt  1  ^  '  X  X 

(o ) 

hence  Z^ '  ^  are  all  normally  distributed  (see  subsections  2  ((5)), 


40 


'2(t6))),  is  distributed  normally 


where 


s.  lO 
la 


=  B  •  X 

.OL  i;a-l,a  ' 
1  ' 


(4.5) 


(4.6) 


=  ?.  •  ''i.S.i’i,.  •  ?«''  ?<. 

1  '  1  1 


(4.7) 


can  be  considered  to  be  the  optimal  estimator  (since,  e.g., 

it  is  both  the  minimal  expected  normed  square  and  posterior  maximum 

{  ci  —  1) 

likelihood  estimator)  of  X .  based  on  data  Z  ,  for  given 

(j) 

partitioning  In  addition,  it  also  can  be  shown  that  „ 

(j)  (j) 

is  the  covariance  matrix  of  error:  ^ia^  •  (See 

Ref.  .125]  for  these  and  other  related  results.  For  convenience,  the 
dependence  on  j,  for  a^j  is  not  indicated  for  the  X,  A  and  related 
terms  below.) 


^i;a-l,a  "  ^ifa-l,a  ^^i  >  Q  > 

M  E(Xj^ 


(4.8) 


=  optimal  estimator  of 

X.  given 
la  ^  1 


and  , 
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,(«*) 


—  Cov(X.  !  Z. 


ia*  X 

=  covttj^, (4.9) 

*=  covariance  matrix  of  error  between 

X .  ,  and  X .  , 

i;a“l»o  1® 

A 

Xi  ^  and  obtained  recursively  as  outputs 

of  the  standard  Kalman  filter  (for  linear  Gauss-Markov  data  measure¬ 
ment  and  target  motion  models  -  which  is  the  case  here;  see,  e.g. 

Ref.  .I26J),  For  completeness,  these  equations  are  presented  below: 


Define  for  all  j  ^  a  ^  0  ,  for  given  Q 


(j) 


S£e(X.  q'^') 

i;o,a  1,®'  1  * 


(®)  n(j) 


and 


df _ (jl  . 

Jttr® 

=  COV(X,  „  1  Z, 


Ai;a,a  “  C°^<V;®,®  “  ^i,®> 


(4.10) 


(4.11) 


1,0  '  1 

Then  in  the  notation  developed  here,  for  given  O^o^j 


A  —  m  ^ 

^i;  0,0+1  *0+1  ^i;o,  o 


(4.12) 


^i;o,o+l  ^0+1  ^i;o,o  ^o+l  ^  ^o+l  ^o+l  ^o+l' 


A  ^ 

i;o,o  ^i;o— 1,0  ^i;o  ^i;o  ^ 


K.  (Q<i>) 
1/0 


A .  ,  •b'^.  ,  if 

i;o-l,o  ^o  1,0  1,0 


^  if  Zi,o=  ^ 


(4.13) 

(4.14) 


(4.15a) 
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.  «  Bi,a  =  ♦ 


>.  =  V,  (Z.  (innovations) 

i;a  1,01 


(4.15b) 


^  z  -  y 
^i,“  i,“ 


(4.16) 


K'^  -R^-^  .B„  +  A.-l  if  2i„  ^  ♦ 

111  1 


^i,*“,«"  Ja 


i,a-l,a 


,  if  2.  =  ♦ 

lot 


(4.17a) 


i;a-l,a 


(4.17b) 


Then  combining  the  results  of  (4.3),  (4.5)  and  (4.16) 


L.„  =  L'  +  lT  , 

la  la  la 


(4.18) 


where 


Bla  =  (q'3>) 


s .  *  log  2  +  log  det  Z . „ 
la  ^  lot 


(4.19) 


and  (noting  s . 

la 


=  non-random  goodness  of  fit  geolocation 
data  term  for  track  set  i  at  t^,  with 
respect  to 

=  s.  (Q^^h 
la 
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(4.20) 


*•1.  ■ 


df 

ia 


I  V 

ia 


la 


=  random  (quadratic  in  innovations) 
goodness  of  fit  geo location  data  term 
for  track  set  i  at  t^^  with  respect 
to  0<a<.j 


It  should  be  noted  that  a  nontrivial  trade-off  exists  between 

the  magnitudes  of  L'  and  L".  .  Thus,  for  the  same  data  set 

la  la 

if  is  a  partition  with  many  relatively  small  (or  one 

point)  track  sets  and  Q  -*  consists  of  fewer  but  bigger  track 
sets,  it  is  not  clear  which  is  larger  in  general:  -2  log  pr ) 
or  -  2  log  pr (Z^ 1 ) .  In  particular,  disregarding  the  inno¬ 
vations  and  the  sizes  of  the  track  sets  (proportional  to 

the  -  which  also  really  play  a  role  in  these  trade-offs 

*  (i ) 

for  Q  ,  the  will  tend  to  be  larger  (in  an  eigenvalue 

**  (i) 

or  matrix  ordering  sense)  than  those  in  Q  ,  due  to  the 

Cov(X^^  I  Z^  Q^^^)'s  being  larger  (less  data  reduction  for 

the  same  target);  thus  the  in  general  may  be  larger  for 

than  for  ,  Yet  this  implies  also  that  the  ^ 

for  will  be  smaller  and  hence  the  L".  's  for  will 

la 

_  ••(i) 

be  smaller  than  for  Q  . 

In  order  to  evaluate  L' .  and  L".  ,  the  determinant  and 

la  la 

multiplicative  matrix  inverse  of  must  be  obtained.  The  evalua¬ 

tions  of  these  quantitites,  as  well  as  the  implementation  of  the 
Kalman  filter  equations  -  (4.14)  -  (4.17),  in  particular  -  are  made 
more  efficient  by  considering  the  relative  sizes  of  s^  ^  and 
m,  for  each  a;  as  a  consequence,  matrix  inversion  and  determinant 
operations  can  be  applied  to  matrices  no  larger  than  min(s^  m) 
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in  dimension 


Specifically,  if  s.  >  m,  then  for  (4.17),  (4.17a)  may 

be  used;  for  (4.15),  (4.15b)  may  be  used  resulting  in  a  simpli¬ 
fication  for  (4.14).  In  addition  identities  for  and 

log  det  can  be  used.  (See,  e.g.,  Ref.  [25],  Chapter  1.) 

Thus,  multiplying  out  the  block  forms  for  the  matrices  involved 
(assuming 


A.  =  (D.  aT^  -  )"^ 

i;a,a  i,a 


\  ,  • (I  +  D.  'A.  ,  ) 

l;a~l,a  m  la  i;a~lfct 


-1 


=  •{!-  (I  +  D.  ‘A.  T  )"^} 

i,a  m  m  la  i;a“l/a 


where 


where 


«  ft 

1,0,  B 


z .  o”  y  • 
i,a,B  ifa,6 


9 


(4 


(4 


(4 


(4 


(innovations  for  target  i  by  sensor  system  3  at  t  ) 


ifeaMia 


n  •  y 

®o,B  *i;a-l,a 


R  •  (I  “  B  -K.  ) 

i“  ®i.a  i“ 

/  such  that 
I  Z.  a'  ^  * 

\h.a.,~  *  ♦ 


(4.25) 


(4.26) 


where 


il=.  «  •  R 


"1  . 


Kl-  (<-27) 


o>-  is  the  Kronecker  delta  function  (i.e.,  6  ^ 

P  f  P  P  r  P 

if  8"  ^  0-j  =  1,  if8'=n. 

Bence ,  (4.20)  becomes 

T  ~  \  '*  V  ♦  A  •  V 

ia  Lu  i;a,6'  iiafQ" 

/l<6',&''<q\ 


=  0, 


I  such  that  \ 

l==i.6'  ’'t'  ) 

K,»B” 


=  V  vT  -R'i-v. 

.  -  x;a,B  otB  i;a,j 


/I4e  <q  \ 

I  such  that  1 

Vi,«.6''y 

(1^8"  <B"  qA 
such  that  \ 

^iaB'  ) 


^aB/'  '’i,a,6" 


(4.28) 
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l"*^  »£i  (A  ) 

fl<  fi%6^<q^ 
I  such  that 

l^iae-*'  ♦' 
Viaa-*'  ♦ 


(4.34) 


where  in  block  form 


la 


^®i;a;B%e"'l 

/l<e'  ,B 


I  such  that 


^iaB 

^iaB 


^  ♦» 
♦ 


with 


(4.35) 


®i;o;6',6-'  ""  ^a'  a" 


A.- 


B'fB"  aB'  aB'  i;a-l,a  oB 


(4.36) 


The  main  difficulties  in  the  evaluations  of  Li  and  Lr 

la  ia 

for  either  relation  between  s.  and  m  are  as  follows; 

la 

For  Sj^^  >m,  eq.  (4.21)^  e.g»»  requires  matrix  inversion  (of 
size  m)  and  eq.  (4.29)  requires  a  computation  of  the  log  determinant 
for  an  m  by  m  matrix  (the  second  term) . 

For  s^^  4  ro,  must  be  obtained,  preferably  in  the  speci¬ 
fied  block  form  of  the  A..  ^  o^'s  (see  eqs.  (4.34)  and  (4.35)) 

X  f  a  ?  p  f  p 

and  log  det  E.  must  also  be  evaluated  (E.  is  of  dimension  s.  ). 

la  la  la 

In  Appendix  A,  two  iterative  techniques  are  presented  which 
can  help  in  resolving  these  computational  difficulties.  The  first 
is  a  procedure  for  obtaining  the  inverse  in  prescribed  block  form 
of  a  positive  definite  matrix  also  given  in  the  same  form  (for  a 
general  number  of  blocks) .  The  second,  similarly,  obtains  the  log 
determinant  of  a  positive  definite  matrix  in  block  form. 
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Some  simplifications  for  the  computations  occur  when 

B  <.  =  B  ^  and  R  „  =  R  i.e.  are  constants  for  all  3  such 
a,e  oO  aB  oO 

that  Z.  „  ♦.  (This  is  not  too  unreasonable  an  assumption,  since 

la  p 

this  still  allows  the  sensor  systems  independent  operations,  but 
with  the  same  general  measurement  characteristics.) : 

Let 


q.  .  y  1 

A  ^  k^q 

/  such  that^ 


V 


1,  jk ' 


Then 


= 

la 


T  '^l 

q.  ‘B  „*R  „*B  „  , 

^lo  aO  aO  oO  ' 


14.37L 


C4,38l 


and  eq,  (4,21)  simplifies  accordingly? 


and  eqs.  (4.24)  -  (4.28)  simplify  slightly.  Also,  (4.38)  and 

log  det  R  =  q.  -log  det  R  -  (4.40) 

^  -a  ^la  a.O 

somewhat  simplify  (4.29).  Eqs.  (4.31)  -  (4.33)  and  (4.36)  also 
simplify  slightly. 

The  special  case  of  one  sensor  system  present  (i.e.,  q  =1) 
shoyld  also  be  noted.  In  this  case,  using  similar  notation  as  be¬ 
fore,  either  ^  depending  on  whether 

^®spectively.  Since  ♦  means  no  geolocation  measure¬ 

ment  is  made  of  target  i  at  t^  by  the  only  sensor  system  present. 
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1 


is  iu>t  calculated  here  and  in  effect 


hi  =  Lr  F  ( 
la  la 


Consider  then  Z.  ,  Thus  s.  -  r 

lax  ja 


la 


For  r.„  >  m  ; 
la 


T  -1 

D .  =  B  ,  R  ,  B  ,  , 

xa  ol  al  al 


with  (4.21)  remaining  formally  the  same; 


X.  =  X.  ,  +  A.  B  ,  R  ,  V .  1 

ija^a  i;a-l,a  IfO^a  ai  al  i;afi  , 


V.  ,=2.,-B,*X.  ,  , 

i;a,l  xal  al  x;a-l,  ' 


i/» 


^i;a;l,l 

“1  -1  T  -1 

=  R  T  -  R  i*B  ,*A.  ‘B  ,‘R  1 
al  al  al  i;a,a  ol  al 


JS'.  =  V  .  ,  *  A .  ,  ,  » A  .  ,  , 

lo  i;a,l  i;a;l,l  l?a,l 


h'.  =  r.  log  2ir  +  log  det  L  , 

la  la  ^  ^  lo 


log  detz .  =  log  det  R  ,  +  log  det  (I  +D.  -A.  .  ) 

xa  ai  m  xa  x;a— i,a 


For  r .  :<  m 
xa= 


F. 

xa 


T 

B,*A.  ,  i»B, 

al  x;a;l,l  al 


with  (4.31)  remaining  formally  the  same; 


X.  =  X.  1  +  A.  ,  •b'^,.A.  ,  'v.  , 

x;a,a  .  x;a-'l,a  x;a'^l,a  al  x;a;l,l  x;a,l 


equation  (4.43)  remains  the  same; 


t4,41l 

C4,42) 

(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 

C4.49). 
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=  A. 

^ia 


=  (R  -  +  B  ,  A  .  ,  B  ; 

al  ol  x;a-l,a  al 


(4.50) 


log  det ‘Z  .  =  log  det  (R.+B*,  A.  ,  Bi) 

^  io  ^  ol  ol  i;o-l,o  ol 


(4.51) 


with  the  equations  for  ^  and  ^  being  formally  the  same 

as  (4.45)  and  (4.46),  respectively. 

The  next  term  considered  in  eq.  (4.1)  ,  -  2  log  pr(ZQ^^^  1  Q^^^)  , 
is  the  goodness  of  fit  of  the  false  alarm  data  to  (or  vice-versa), 

0  ^  i  j 

Note  first  the  relations  (a  consequence  of  the  assumptions  in 
subsections  2  ((3))  and  2  ((6)): 


=  n  n  . 

0<a<j  \ 

I  such  that! 
\f  Oak  *  ♦  / 


(4.52) 


where  if  Z-  ,  <ji,  f  ,  >  0  ,  then 
Ook  '  ^ '  ak 


(2Tr)  (det  M^j^) 


(4,53). 


^Ootk  —  ^  ^^0aku“®ak^  ’^ak  ^^OakuT  ®oik^ 


^"Ook  ^  "ook  ' 


(4.54) 
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^'^Ook  ^  ^'’Oak  ^^Ook' 


ok 

~  E  (Zoak  -*0ak>’’-^i<"0ak 
(0  =  1 


“  ^Oak^ 


=  ^-^'^Oak'^i’  < 


(4.55) 


^Oak  —  £  ^^Oak  "  ^Ook^  ^^Oak  "  ^Oak^  ' 
(0=1 


■ak 


(4.56) 


'Oak  —  ^^'^^ak^  ’  2  "Oakco 


'ak 


(4.57) 


(0=1 


^■^Oak  ^"Oak  ^^Ook' 


(4.58) 


df  «p  -1  - 

^ak’^^Oak"  ®ak^  ”ak  ^^Oak  "  ®ak^ 


(Equation  (4.54)  follows  from  the  standard  procedure  of  adding 
and  subtracting  within  the  sum.) 

^Oak  J^®Pi^esents  the  random  (geolocation)  fit  of  the  false 
alarm  set  for  sensor  system  k  at  t^  with  respect  to  , 

0  <  a  <  j 

L'rt  1,  measures  the  scatter  of  the  false  alarm  data  for  sensor 
Oak 

system  k  at  t^  with  respect  to  Q^?  '^Oak  unnormalized 

sample  covariance  of  the  false  alarm  data  and  is  the  sample 

mean,  for  sensor  k  at  t  . 

L^Q^k  measures  the  bias  between  the  observed  and  predicted 
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^Oak  ”'®^sures  the  total  goodness  of  fit  of  the  false  alarm  set 

for  sensor  k  at  t^  with  respect  to 

Computational  problems  for  the  false  alarm  scores  appear 

minimal;  the  real  difficulty  lies  in  the  modeling  -  specifically 

in  the  choice  of  0  's  and  m 

qK  aK 
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The  third  teinn  considered  in  eq.  (4.1)  -  2  log  pr  (Y )  , 

is  the  likelihood  (or  goodness-of-fit)  of  the  non-geolocation  target 
attribute  data  with  respect  to  . 

It  follows  from  the  assumptions  made  in  subsection  2  ((7))  that 


(j) 


n 


I  such  that! 

\’'iak  *  ♦  j 

Then 

-  2  log  pr(Y^^^  )  , 

=  -21og  5^  pr(Y^^MH^^\  •pr(H^3)) 

(over  all\ 
outcomes  1 

of 


I  such  that  I 

\y,<3> 


(j) 


where 


(4.62) 


(4.63) 


—  -  2  log  pr(Yj,  |Q^^^  ) 

=  -2  1og(  53  pr(Y.  ^^^1  H.  'pr(H.)) 

s  4  «  V  X  X  X 


(over  all  \ 
outcomes  of  I 
H,.  C  ) 


(4.64) 


In  addition,  note  that 


pr( 


\y(3>/  Aiii”  \ 

^  »  I  4-V»a4‘  I 


(4.65) 


such  that 

Z.  ^ 
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(4.66) 


pr(2.  = 


0^o^j\ 
such  that) 

^  ^  J 

la  y 


The  goodness  of  fit  of  overall  data  density  at  time  t^ 
with  respect  to  is  analogous  in  form  to  that  of  the  geo¬ 

location  Only  data  given  in  (4.2)  et  passim,  and  can  be  decomposed 
into  a  product  of  two  factors,  one  representing  geolocations  and 
the  other  non-geolocation  attributes: 

pr ( Zi J  Zj ,  q'Z>) 

=  (4.67) 


(provided  Z .  y  4>  and  Y .  y  <f>  for  some  0  ^  a  ^  j ) 


Also 


(j)  -  - 


where 


2  log^  n 


pr(Y^^|Y.  Q^^^) 


70  <  a  <  j  \ 

I  such  that! 

\y,  7*4  y 

Z)  .-2  log  prCY  |Yy“-^>  ,  Q*3>) 

/0<a<j  \  ^ 

f  such  that! 

/o<?<  j  \ 

I  such  that! 

\Yi 


=  L  .  (Y. 

la  1  '  ^  ’ 


=  -2  log  pr  (Y^^|y^^“”^\ 


(4,68) 


(4.69) 


noting  for  0  ^  a  ^  j 


i 


pr(Y.  -  .  £  .pr(Y.  1H.,Y  ^®“^^Q^^^.pr(H.) 

^  .  /over  all\  i  i  i 

f  outcomes  \ 


\of  n^ecj 


(4.70) 


If  a  and  b  are  relatively  small  integers,  by  using  the 
(known)  distribution  functions  of  the  I/.  ..  's  or,  equivalently, 

i  jK 

the  p(y"’,Y-^;Q^^^ ) 's,  and  eq.  (2,4)  (see  subsection  2((7)),  then 
equation  (4.63)  can  be  evaluated.  However,  for  relatively  large  a 
and  b  this  may  require  -  because  of  the  discrete  nature  of  the 
distributions  involved  -  a  large  number  of  operations.  Consequently 
the  following  approximation  is  proposed  for  the  modeling  of  the 
random  variables  either  relatively  large  a 

and  b,  and/or  when  the  non-geolocation  attribute  set  C  is  per¬ 
haps  better  modeled  as  a  contiguous  subset  of  b-dimensional  Euclidean 
space : 

In  equation  (2.3b),  assume  each  ^ij)^  (for  /  (j)  is 

normally  distributed  as  N  (0,  where  is  the 

I3K  IjK 

submatrix  of  R,  corresponding  to  outcome  T. ..  with  respect  to 
{1,2,..,  b},  R  being  a  fixed  (positive  definite)  matrix  of  dimen¬ 
sion  b,  to  be  determined. 

Define,  analogous  to  Rq  , 

i 


^iaq 


I  For  only 
1  ^  k  ^  q 
such  that 
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Replace  equation  (2.4)  by: 

is  normally  distributed  as  Nj^(E(H^),  Cov(H^)),  with 
E(Hj^)  and  Cov(Hj^)  to  be  determined  in  relation  to  attribute  set 
C.  The  following  trivial  identity  is  used  here: 


H.  t  \  =  H.,  for  all  times  t  , 

X,  id)  X'  a 


(4.71) 


From  now  on  for^^iven  time  t.  '  and  ,  the  function 

dependency  of  quantities  on  will  be  often  omitted,  unless 

cunbiguity  results  or  emphasis  is  desired. 

Then  to  evaluate  a  Kalman  filter  can  be  applied  to 

(2.3b)  and  (4.71)  yielding,  analogous  to  the  computations  for 
L.  (see  egs.  (4.5)  -  (4.20),  for  0  4  ct  £ j »  for  given  : 


A  A 

H.  =  H. 

i;a,a+l  xjUfa 


(4.72) 


^i;o,  a+1 


^i;a ,a 


(4.73) 


H.  =  H.  ,  +  K-  ’v. 

x;a,a  x;a-l,a  x;a  xjo 


(4.74) 


T 

* 


x;o 


L.  ,  ’B  ’Z.  fif  Y.  A 
i;a-l,a  .a  l,a  '  Xa  ^ 


,  if  Y.  =  (j. 

xo  ^ 


(4.75a) 


i'*  T  -1 

A.  •  8  .  R 

I  i;a»a 


,  if  . 


,  if  Y.  =4 
la  ^ 


(4,75b) 


a 


(8^*R  ^*8  +Ar^  ,  )"^  ,  if  Y.  ^ 

x;o-l,a  Xa  ^ 


^i;a-l,o 


,  if  Y.  =  (j. 
xa  ^ 


(4.76a) 
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i;a-l,a 


,  if  ^  ♦ 


'  ’'ia  “  ♦ 


(4.76b) 


2£  y.  -8  .li.  , 

i?o  ^  io  .o  i;a-l,( 


d'f  ** 

Z.  .  B  -A 


8„  +  R„  . 


W  ••  •  ^  S»  ^ 

xta-l,a 


Then 


^ia  =  iia  +  ^ia 


where 


df 


i " .  0  -A .  •  log  2ir  +  log  det  E 


la  xa 


ia  * 


(4.77) 


(4.78) 


(4.79) 


(4.80) 


(4.81) 


(4.82) 


Analagous  to  the  cases  ®i<x  ^  ®ia  evaluate  the 

above  equations,  we  roust  consider  separately  the  cases  ^  ^ 

and  A  £  b  : 

For  A>  >  b,  use  eqs.  (4.76a)  and  (4.75b),  obtaining 

Xd 

1  —1 

A.  *  (V-  +  A.  ,  ) 

x;o,a  xa  x;a-l,o 


"  ^i;o-l,a  ^^b  ^  ^io’^i?a-l,a^ ' 


(4.83) 


etc. , 
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9 


(4.S4) 


Pj.  SI  8  T.R-1b 


/  S  V  iae^ioB^iae 
/I46 < q  \ 

at\ 


I  such  that^ 

V'^iuB  ^  ♦  / 


A 

H 


i  *a»a 


^i;a  ,a 


A^6<q  \ 

I  such  that  A 


(8?  ) 

laB  laB  i;a,B^ 


\’'iaS  »'  ♦  J 


(4.85) 


&V  R  6 

i;aB  ^  iaB  ioB*” 


laB  XaB  i;a“l,a» 


*0  •  6^*  8"^ 

'  '  '  /l4S",r"<q 

/  such  that 

(  ^iaB'  ^ 


A. 


df 


i;a;B^B-'—  r ‘’^iaV-^i”aV  *®iaB  ' '^i;aa '^iaB^ ’’^i'aV' 


L"' 

la 


T  ^ 

'  2^  ,''i;aB'^i;o,B%B^ ''i;aB 

such  that  I 


Y._„.  4>  / 


ioB- 

r 

iaB‘ 


f  ^  j 


(4.86) 

(4.87) 


etc. 
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for  1  ^  Y  i  b 


for 


where 


where 


for 


and 


Hence , 


f 


/ 

/  such  that  Y.  „ 

/  Ia$ 

I  and  there  is  an  n, 

I  1 ^  n  ^  b,  such  that 

V  ^  ^ioPn  ^  ^iaP 


liY  <  b. 


(4.94) 


(4.95) 


(4.96) 


df 


'iaPY 


,  If  r  / 


(4.97) 


(4.99) 
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Also^ letting  in  scalar  form 

=  (A.  -  4 

i;o,a  '  i;0/“;Y  „ 

*  i  T  /  T 


<  b 


(4.101) 


then 


®  •  Of  «  ^  o,  o,  *  ®  •  «  fl  -  /T 

lOp  !;<*,«  rag  =  (A 


i;a,o;Y  ,Y 


’'''Vb- 


Note  that 


iae 


log  det  =2  '"b.  . 

n=l  io6;n 


=  2  log 


YeT.  o 
ia0 


For 


A.  =  (I,  -A.  ,  ‘F.  )  'A.  , 

x;a,a  b  i;a-l,a  la  i?a-l,a 


where 


df  2 

(liB',B~  <q,\ 
such  that  \ 


ia 


^  v^iae'^ioe'^g"  ®iog 


(4.102) 


(4.103) 


(4.104) 


(4.105) 
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I . 

i;a,a 


L  ■  1  •  >  ^ 


1  ^  6'#  B"  <  q 
such  that 

’'io8'  ''  ♦  ' 


such  that 
^iaB'’  ^  ♦  ' 
^ioB^  ^  * 


(4.106) 


(4.107) 


df  T 

®i;a;B/,B'''  ~  ^B",  B"'’*^iaB"  ‘‘’  ®ioB'**ifa-l»a’®io 


(4.108) 


Further  simplifications  can  be  carried  out,  since  R  is  in 


diagonal  form. 


For  example,  if 


^i;a;  B^f  B"  “  ^^i;a;B'’,B'^;n'#n^ 


li 

li 


(4,109) 


then 


®ioB'’^iaB'B"*®iaB'"  “  ^^i,a ,  B B "7 y' » Y ^ 


1  <  Y'.Y"<b 


(4.110) 


where 


^ioB',B"^Y^y 


if  ^  "fiaB'  'fiaB 


^i;a;B"f  B^,n',n-'  ' 


^  ^ioB'n'  ^  ^iaB^  ^ 
Y"  =  ^idB-'n"^  '^iaB"  ' 


(4.111) 
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SlnlXarly, 


a; 


laB 


-•*l,«,a%6~'’i,aB“np  =  ) 

VioB'r'b/ 


(4.112) 


where 


^ioB'B^Y  " 


0  ,  if 


ioB 

n%  n"  n 

r\^  =  1 


,  if 


^  =  Ke'n'  "  ^i«S 


(4.113) 


Alternatively,  computed  directly  using  eq.  (4.71) 

and  the  fact  that  is  normally  distributed  with 


E(y.  (j)lQ(j))  =  8(j) 

•E(H.) 


(4.114) 


and 


.xT 


Cov(Y.  •Cov(H.)  *8^^^ 


(4.115) 


where 


(j)  df 


/h 


B. 

O 


(4.116) 


Let 


(j)  df 


E  4, 


a=0 


la 


(4.117) 
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and  assxime  Then 


i.  (j)  ^  L  '  +  1  tj)- 
i  1  ^ 

^3^  2£  ‘log  2ir  +  log  det  Cov(Yj^ ) 


(4.118) 


(4.119) 


^^(j)  df  (Y^Cj)  _  |Q^3)  j)T^^^^^Y^(j)  |Q(j)  j  ^Y^(j)_g^^^(j)  |^(j)  ^  j 


..  -1 


(4.120) 


Analogous  to  eq.  (4.29), 


logdetCov(Y^^^^  |Q^^^)  =  log  det  +  log  det  (I  •Cov(H.)),  (4.121) 

4  ^  ^ 


pO).  df  8(.jl’'.pU).-\Btj) 


=  t  "i, 


(4.122) 


CoylXi  =  (B^^^  Cov(H.)  8^^^  +  R^^^) 

'■■■■*  i  1  • 


f.-l 


_1  _2^  qt  _ 

8^^^  •  (Cov(H.)  +P^^^) 


1  1 


(4.123) 


where 


jj(-j)  ^ 


(4.124) 
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-E  (Y 

i  i  i  i  ^ 


svK;h  "that 

^iak  ^  ♦  C4a25) 


E(Y^^1q^^^)  =  8^«EtH^)  , 

and 

E(Y.„^|Q<1'  =  ,-  l<k<q 


(4.126) 


t4.127) 


Then,  can  be  evaluated  by  using  (4,122)  -  (4,127)  in 

(4.120) . 

If,  furthermore,  R  and  Cov(H^)  are  diagonal,  additional 
simplifications  can  be  achieved  by  use  of  (4,90)  -  (4.95),  replacing 

~iia,6  ‘’y  ’"ic.r  (■'•W  -  (4.100). 

An  approximation  for  ,  the  computations  of  which  do  not 

depend  on  E(Hj^),  which  is  more  accurate  for  larger  (in  the  positive 

definite  ordering  sense)  Cov(H^),  is  given  in  Appendix  B,  Thus, 

using  the  approximation  no  knowledge  by  the  observer  is  needed  of 

E(H^)  and  Cov(Hj^)  ,  except  that  for  all  i,  Cov(H^)  exceeds  a  certain 

fixed  large  bound. 

\^en  Cov(Hj^)  is  diagonal. 
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5,  SOIME  DISTRIBUTIONAL  .  PROPERTIES  OF  THE  SCORE 


As  in  the  last  subsection,  it  is  assumed  that  the  prior  dis'^ 
tribution  of  is  uniform  so  that  J-  and  I""  need  only  be 

considered. 

*  Then  from  eqs.  (4.1),  C4. 2),  (4.4),  (4.18),  (4.19),  14.20), 
(4.60),  C4.61),  (4.59).  (4.54)  and  (4,62)  -  (4.64)  and  the  assumptions 
made  in  subsection  2  (J'' (Q^^^  ,Z  )  ,  as  a  random  quantity 

functionally  dependent  on  is  distributed  as  the  sum  of  a 

constant,  a  chi-square  random  variable  and  a  statistically  independent 
discrete  valued  random  variable.  Specifically, 

J- (Q^^^  ,Z  ^^^ )  =  »  Z^^^)  +  (5,1) 


is  the  constant term; 


a<i<M<3) 
such  that 


Oak 


0  4  a  ^  j 


^1  ^  k  ^  q 
such  that! 

^'Oak  ^ 


(5,2) 


(5.3) 
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(since  all  (Lr  1q^^^)'s  are  distributed  as  Xe  statistically 
1  s . 

la 

independent  for  all  (i»a) *s  and  statistically  independent  of 

all  ) 's  which  are  distributed  as  )  is  dis^ 

.  2  ,  ^o»k  ^ak 

tributed  as  Xp  n  where 

j 


(5,41 


(5.51 


has  a  discrete  distribution. 

Under  the  normal  approximation  made  for  (see  the 

results  follov/ing  eq.  (4.70),  especially  (4.113)  -  (4, 115)), it  follows 
that  the  discrete  valued  random  variable  ,  Y^^^)  is  replaced 

in  the  sum  comprising  ^^^)  (eq,  (5,1))  by 


/I  <  i  <  M 
I  such  that 


(j) 


+  L 


(j) 


V 


(j) 


(see  eq.  (4.113)), 
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which  causes  an  adjusted  computation  for  J'.  J' ,Z ) 
is  now  distributed  as  the  sum  of  a  constant  and  a  chi-square  random 
variable.  Specifically, 


j-  ,Z  ) 


Jf  , 


where 


is  the  constant  term  ; 


(5.6) 


(5.7) 


(Q 


(j) 


df 


(5.8) 


(5.9) 


It  should  be  noted  that  if  0^3^  the  distribution 

vu  ly  ;|Q  -  o  )  IS  in  general  not  obtainable  in  a 

simple  closed  form  -  even  under  the  normal  approximation  made  for 
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This  is  due  mainly  to  the  mismatching  of  covariance 

matrices  of  formally  assumed  data  model  by  the  correlator  tracker 

through  Q  ,  with  respect  to  the  actual  data  model  through  Q  , 

For  example,  in  reality  let  q=l^  j  =  5, 

P  hO,  B=B,=  B.,  ,  R  =  R  -  =  R-,  .  (See  eqs.  (2,1)  and 
o  U  o  al  01  al  01 

(2.2).)  Thus  no  false  alarms  nor  nonr-geolocation  attribute  data 
are  truly  present  and  all  observations  and  target  motion  are  homo- 
geneous  and  stationary. 

Suppose  also  that  only  (geolocation)  data  Z,  , 
for  a  =  0,  1,  5  =  j  is  observed  (see  the  concluding  part  of 

subsection  2  ((8))  for  use  of  notation),  where  for  each  a,  in  general 
“  ^ial'  ^  ~  unknown  positive  integer. 

Suppose  now  the  observer  uses  a  correlator- tracker  determined 

® (5)  *  f 

by  partitioning  Q'  \  consisting  of  a  single  component  Qj  , 

*  (5) 

=  {  (l,a,l)  I  a=  0,  1,  ..,5},  corresponding  to 

Z.  =  {Z,  J  a  =  0,  1,  ..,  5}, 

±  X ,  a  /X 

/ c>  0  / C\ 

Now  Q'  ^  =  Q'  '  really  holds,  then  without  loss  of 
generality,  Z^  ,  =  ^lal'  ”  ~  observation  and 

target  models  in  eqs.  (2.1),  (2.2)  combine  to: 


^101  \ 

/BOI  \ 

1  ''loi  \ 

"ill ' 

'  ®01-*0  ' 

^101 

''ill 

• 

"l5li 

• 

1 

i  ''151  / 

■  —V  —  ■ 

(5.10) 
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yielding  being  normally  distributed) 

Cov(Z^^Mq^^^)  =  ^•CovCXj^Qj^)*^’’  + 

Now,  has  random  term  {see 

eq.  (4.20)). 


L3  =  L5 


(5.11) 


a=t 


''12  •^12' 12 


=  (Z^^^  -  E(Z^^^)'^(Cov(Z^^^“^  ‘(Z^^^  -E(Z^^h) 


*  (5) 

(conditioning  here  is  on  Q  ) 


Then  (5.10)  and  (5.11)  imply  that  (L3(Q^^\z^^^ )  ) 

2 

distributed  as  Xr  »  where  now 


IS 


«5=  E 


s . 

la 


(5.12) 


#51  00  (5) 

On  the  other  hand,  suppose  Q  =  Q  really  holds,  where 

now 

=  {  (l,a,l)|  a  =  0,  1,2}, 

=  {(2,a,l)  I  a  =  3,  4,  5}  . 
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Thus  corresponds  to  =  IZj^^j^|a=0,l,2} 

00  /  C  \ 

and  Qj'  *  corresponds  to  ^^20!^®“  3,4,5}  =  { j a=3 , 4 , 5} . 

This  yields  from  eqs.  (2.1),  (2.2)  the  combined  observation 

model: 


But,  since  here 


e(z<s>|q'=>) 

\6-.E(X23i)) 


,  and  most  importantly 


Cov(Z  )  = 


|^'Cov(X233^)^' 


+  (5.14) 


#  c\  *  f 

does  not  match  in  shape  (nor  size)  Cov(Z'  iQ-  M  which  is  the  central 

factor  in  the  quadratic  fom  in  (5.11),  it  follows  that 

(Lg (Q , Z I  Q^^^)  is  not  even  distributed  as  a  noncentral  chi- 

0(5)  (5)  00  ( 5) 

square  random  variable;  hence  neither  can  (J'(Q  ,Z  ) [Q  )  have 
a  simple  distribution. 

Note  that,  by  similar  reasoning,  if  a  correlator-tracker  is 

»» (5)  "  (5^  (5}  I  (5) 

used  which  formally  assumes  Q  to  be  true,  then  (L(.(Q  ,Z  )lQ  ) 
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is  distributed  as  x?  «®  in  (5.12),  but,  as  before,  (L, 

has  a  complicated  (non-chi- square,  non-noncentral  chi-<-square) 
distribution  in  general. 

Consequently,  in  general,  for  Q'-*'  ^  Q  ,  since  we  cannot 
obtain  a  simple  computable  distribution  for  J' (Q , Z ) . 

similar  remarks  hold  for  the  difference  ^  J' (Q^^^  , Z  J' ,Z  ^ )  f 

conditioned  on  Q  ,  Thus,  the  computation  of  a  threshold  T  Csee 
eqs.  (3.8)  -  (3.10)  based  only  on  a  given  significance  level  of 
discrimination  B  (obtained  by  solving  B  =  Pr  (Decide  ^^2) 1^(1)  ~ 

Pr(Jj>T|Q  for  T)  appears  equally  infeasible. 

* (i)  * (i) 

In  summary,  J'(Q  •'  ,Z  "^  )  can  be  used  two  ways  in  a  real- 
world  situation; 

(1)  For  determining  how  well  a  given  correlator^tracker , 

through  its  partitioning  of  data  really  fits  data 

by  evaluating  the  cumulative  distribution  function of  the  random 
quantity  (J  (Q^^\z  )  at  the  outcome  point  J  (Q^^^  ,Z^^^  ). 

This  distribution  in  general  should  be  computable,  since  the  random 
quantity  here  is  the  statistically  independent  sum  of  a  computable 
constant,  a  chi-square  random  variable  with  a  computable  number  of 
degrees  of  freedom  and  a  discrete  random  variable  which  has  computable 
characteristics  (theoretically  at  least,  if  parameters  a  and  b 
are  small)  (see  eqs.  (5.2)  (5.5).)  (An  approximation  based  on 

simplifying  the  computational  burden  for  the  discrete  random  variable 
term  is  given  in  eqs.  (5.6)  -  (5.9). 

(2)  For  comparing  two  (or  more)  given  correlator-tracker  schemes, 
operating  on  the  same  data.  In  general,  up  to  an  adjustment  for 


74 


different  decision  losses  or  gains  and  presence  of  prior  distributional 
information  on  (see  eqs.  (3.6)  -  (3.10)),  that  correlator 

tracker,  determined  by  partitioning  ,  is  chosen  among  a  given 

set  say,  for  which  J '  (Q  ^  ,  Z  ^ )  is 

minimal.  However,  the  statistical  significance  of  the  differences 
between  the  values  of  for  different  correlator- 

trackers  is  apparently  difficult  to  determine. 

Note  that  using  the  distributional  results  in  eqs.  (5.2)  - 

(5.5) ,  for  example,  I‘'(Q^^^,  is  easily  obtained  as 

q<3))  =  (Q(j))  +  Sj  +  E(J3(Q^^^  Y^^M).  (5.15) 

If  the  normal  approximation  is  made  (for  lo^^^)  ,  then  eqs. 

(5.6)  -  (5.9)  imply 

=  J£'(Q  (5.16) 

i'  evaluated  as  in  (5.16)  can  be  used  also  as  measure  of 
average  performance  of  Q  (with  respect  to  averaging  the  data 

Z^^^)  and  can  be  used  analagously  to  J'  in  (2.)  for  comparing 
average  goodness  of  fit  of  several  correlator-trackers  in  question, 
by  choosing  that  one  minimizing  the  corresponding  value  of  i' . 
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SUMMARY 

A  simple  scoring  rule  for  correlator- trackers  is  developed 
in  this  report.  Mathematical-logical  justifications  for  the  use 
of  this  score  are  demonstrated.  Detailed  computations  necessary 
to  implement  the  rule  are  exhibited,  along  with  a  determination 
of  its  statistical  distribution,  a  form  of  the  chi-square. 

In  the  second  part  of  the  study,  numerical  examples  will 
be  given  illustrating  the  suitability  of  the  scoring  technique 
for  use  in  surveillance  in  a  real  world  setting  and  leading  to 
sensitivity  analysis  with  respect  to  the  key  parameters  involved 
in  correlation.  Future  work  will  concentrate  on  both  extending 
the  applicability  of  the  score  and  obtaining  further  analytic 
properties . 
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Appendix  A 

Matrix  Inverse  and  Log  Determinant 
of  a  Positive  Definite  Matrix  in 


Prescribed  Block  Form 


Let  G  =  (G..)  be  a  given  positive  definite  matrix  in 

1  i  j 

block  form  with  all  G^^  assumed  square.  Thus  all  G^^'s  are  posi- 

T 

tive  definite  with  G..  =  G.  .  ,  etc. 

13  3,1 

Let  G  ^  ^  (A..)  ,  with  the  A..'s  to  be  determined. 

1  <  i.j  .<  n 

Then 


1.  For  1  £  i  ^  3  £n 


k=i  +  l 

assuming  the  sum  is  zero  for  i=n,  where  for  any  i£i,3,  k^n, 

G  0)  „  G 
^ik  ^ki 

For  i  =  1 ,  compute  for  j  =  l,2,  j^k^n. 


(A.l) 


C  (1)  df  . 

=  Ckj 


(A.2) 


For  i  such  that  n^i>l,  compute  for  j  =  l,i+l,  j£  k£  n. 
(For  i  =  n,  let  k=j  =  i  =  n,) 


G  -  G  .  G.(i-P 

kj  jk  .  k,i-l  1-1, 1-1  1-1,3 


(A.  3) 


A-1 


2. 


log  det  G 


(A. 4) 


n 

=  53  log  det 

i=l 


)  . 


For  i 


such  that  n  >  i  >  1 


compute 


for 


i  +  1  ^  j  i  i  • 


(For  i®n,  let  k=j  =  i  =  n.) 

The  proofs  of  the  above  results  easily  follow  by  successively 
applying  the  matrix  and  determinantal  identities  in  Ref.  [25], 
Chapter  1,  pp.  32,  33. 


A-2 


Appendix  B 
An  Approximation  to 


Assume  rank  Define 

i 

1  111* 


where 


Let 


^  i  i  i  i  i  i 

~  mineig  (Cov(H^)) 

^ijl  ^  inineig  Cov(Rj) 

^ij2  ^  mineig 


Then,  for  any  arbitrary  given  e  >  0,  if 

hi  >  (l/h.j2)-max(0,(  ||  Y .  ^ ^ ^ -E(Y.  ^  |1  Ve-h.^j)-l), 


then  (see  eq.  (4.120)) 


0  < 


<  e 


Proof : 

Note  that  by  using  matrix  identities  (Ref.  [25],  Chapter  1;  see 
also  eq.  (4.123)) 


B-1 


where  (see  eq.  (4.117) 


T 

fjCi)  =  g(j),(-p(j).Cov(H.)*P^^^  .  (g.2) 


(B-I) 


Then 


Ra)-\n(j).R(j)-i 


=  .(R 

i  i 


1  1 


)*R 

i  i  i 


(j) 


Since 


i  i  i 


•  -h 

^  maxeig  (R^^^  ^.j 

i  i  i  *  P  ^ 


(B.3) 


(B-4) 


(in  positive  definite  matrix  ordering  sense) . 
Now 

maxeig  (K  n  r  -» 

i  i  1  ^ 

-1 


11  1 
-1 
-n^ 


=  maxeig  (P  ^  ^  Cov  (H .  )  .p^j^) 

i  ^  i  ■  ' 

=  maxeig  (Cov(iL)*  +  l^) 

=  1/Cmineig  (Cov(H. ) -P^^^)  +  1) 

1  • 

1 

<  l/(hi.h..2  ^  1)  , 


(B-5) 


B-2 


Noting 


that  =  0,  it  follows  that 


(j) .E(Y. =  H. -8^^^ -ECH.)  =  0  and  thus  combining  (B-1)  - 

1  1  1  i 


(B-5) 


0  <(Y.(S’-E(Y/Jb)'';(Cov-l(Yj«b-Hi«’)-(Yifi’-E(Y/ih 


(omitting  the  conditioning  on 


,1.0)".  i.O) 
1  1 


~h  ■  ~h 

<  ||Y.^^^  -E(Y.^^^I1  ^‘maxeig  *  (1/ (h.h.  .  ,  +  1))  •  I ’R^^^  ) 

1  1  j  ^ 

A 

<  IIYjH’  -  E(Y.O)||  -a/fhi.iChjhj.j*!))}  .  tB-6) 


B-5 


Appendix  C 

Flow  Charts  for  Computations  of  the  Score  J' 

The  number  of  sensor  systems  is  q,  known  and  fixed. 

At  each  time  t.,  new  raw  data  (possibly  vacuous)  Z.  is  ob- 

3,  VP  >  j  3 

served  before  is  determined.  Z.  =  (Z  .  ^) 


'y  ,3  ,k' 


Y  1  y  .  a  y  TH  ■ 


jk 


^  1 


For  each  sensor  system  k,  at  t,,  m.,  ^0  is  the  known  number  of 

data  reports  2  .  ,  observed.  (When  m.,  =0,  data  becomes 

Y  » 3  » ^  3  ^ 

missing,  i.e.,  Z  =  <(1  when  m. ,  >.1,  only  at  most  one  report  can 

Y3  K  3K 

correspond  to  each  true  target  -  the  remaining  are  false  alarms.) 

Also,  ^yjk  decomposed  into  a  possible  geolocation  (target  or 

false  alarm  -  not  known)  data  component  which  if  non- 

vacuous  is  r^j^  by  1,  r^j^^l  known,  and  a  possible  non-geolocation 

target  data  component,  Z". ,  ,  which  if  nonvacuous ,  rules  out  the 

y3  ^ 

associated  geolocation  component  being  a  false  alarm,  and  is  of 
dimension  ^b;  C  is  the  fixed  known  non-geolocation  target  attri¬ 
bute  set  of  b  by  1  vectors,  each  representing  a  feasible  evaluation 
of  b  given  attributes. 

For  each  sensor  system  k,  at  time  t.,  B.,  is  known  r., 

3  3k  3K 

by  m  geolocation  target  measurement  matrix,  and  R.,  is  a  known 

3  ^ 

corresponding  r^j^  by  r^j^  positive  definite  measurement  error  co- 
variance  matrix.  m  =  dim  (X.)  is  known;  X.  is  unknown  i^^ 
target  state  parameter  vector.  See  Eqsi  (2.1)  and  (2.2)  for  further 
explanations  and  other  related  definitions. 

Prior  to  t ^  ,  partitioning  and  target  state  parameter 


C-1 


vectors,  one-step  predictions  from  the  Kalman  filter  are  available; 

X.  .  ,  state  estimator;  and  A  ,  covariance  matrix  of  exti- 

mator  error;  i“  1,2,.. 

Also,  available  prior  to  t^  are  geolocation  scores 
i  =  0,1,,.  (i*0,  corresponding  to  false  alarm  set;  i^l  corres- 

ith  _ ^  ^ _ _ _i_ _ ^  _ * _ j  . _ 


ponding  to  i  target  track  set,  determined  by 
.attribute  scores 


Q''-’-'),  non-geolocation 

^  ,  i  =  1,2,...,  possible  Kalman  filter  pre¬ 

dictions  for  ^  true  b  by  1  attribute  vector  of  target  i  -  as 

A  A 

1  T*  1  -1^  etc. 


Also  available  prior  to  t^  are  the  overall  geolocation  target 


data  score,  -2  log  prCZ^^^  overall  false  alarm 

(geolocation)  data  score,  -2  log  pr  (Z^  (O  and  the 

overall  non-geolocation  target  attribute  score  -2  log  pr ”^^) 


Following,  the  reception  of  nevr  data  2.  ,  based  on  all  of  the 


old  data  Z  and  the  new,  combined  into 


and  based  on 

possibly  old  partitioning  new  partitioning  of 

is  determined  by  the  observer. 

Once  is  determined,  then  the  total  data  can  be 

decomposed  as  =  {Q^^^^li  =  0,  1,2,,..},  where  the  i^^  track 

set  (target,  for  i^l,  false  alarms  for  i  =  0)  is 

=  ^i^^^  ~  l£k<q}.  The  are  the  same 


as  the  Z^^j^'s  with  the  i’s  and  y's  replacing  each  other,  except 


^Oak  “  ^Oqk  "  ^^Oako) 


l^u^fjj^}^  the  false  alarm  set,  where 


fj^^O  known  relative  to  0^^^ 


(All 


Oaktu 


’s,  as  well  as  all 


Zink's,  unless  vacuous  are  ^  r^^^  by  1.) 


C-2 


=  ,  in  which  case. 


Some  of  the  ^iak'®  vacuous 

those  are  generally  marked  so,  deleted  from  the  set,  and  do  not 
contribute  to  the  computations  (actually  contribute  zero  in  value 
to  the  various  sums  calculated).  Relative  to  i  is  formally 


assumed  known  and  Z.  , 

lak 


’'iak 


‘iak 


is  the  r^j^  by  1  geo¬ 


location  data  component  and  Y  is  the  non-geolocation  target 

attribute  data  component  -  a..,  by  1,  Furthermore,  random  set 

IJ  K 

^iik  ^  ^^iik  1’  ^iik-a  }c{l,  b}  is  observed,  and  may 
be  vacuous,  in  which  case  a..,  =  0.  Without  loss  of  generality, 

1 1  K 


l<-b..,  -•^..<b..,„  ^b.  Associated  with  the  Y.  ,  ’s  are  measure- 

-  ijk,l  ijka^jj^-  ictk 

ment  matrix  8.  ,  and  error  covariance  matrix  R.  ,,  (see  remarks 
lak  lak  ^ 

between  eqs,  (4.70)  and  (4.71),  and  eq.  (2.3)  and  following  discussion.) 


Note  also  the  notation:  Z 


0)  =  z^Cj)  ,, 


(Z, 


alarm  set  up  to  t.;  Z, 

J  ’  + 


Cj) 


0  ■  "0  ^“Oak 

(.i) 


0  ^  ot  ^  j  };  false 


=  {Zj'"*'^  1  i^l},  for  the  set  of  all 


geolocation  data  Z^^^^  =  t  ^  ^ak  ^  ^  ^  ®  =  ^^iak  track 

set  i;  up  to  t ^ ,  for  all  i^l,  and  Y^^^  =  ^1),  for 

the  set  of  all  non-geolocation  data  of  track  set  i,  Y^^^^  = 


l^h^q;  O^a^j  f  <l>)},  up  to  t  ^ ,  for  all  i>l. 

For  additional  clarifications  and  definitions,  see  the  main  text, 
especially  subseet-i-on  2  of  the  Analysis  Section. 

The. gei}e.ral:  convention  assumed  tin  these  flow  charts  is  that  any  sum 
whose,  index  ,  set  of  jsummatidn  is  the  empty  set  ((^),  is  set  equal  to  zero. 


C-3 


Define  as  the  first  j  >^0  at 

sampling  time  tj ,  for  which 


^  ji<  ^  .  Hence, 


i  ♦ 


C-5 


Geolocation  Target  Data 


C  ~  {ill<i  6  f  (t>  S  i  /  C..UC.  UC.._} 

=  set  of  all  track  sets  consisting  of  two  or 
more  points  which  are  obtained  by  breaking  up 
previous  adjudged  false  alarms,  given  ^  ^ 

This  is  equivalent  to:  ^  ?  , 

i.e.  are  some  previous  track  sets  broken  up?  _ 


-2  log  pr 

=  -2  log  pr(Z^^j'^^ 

+  Y'h.. 

CieC.jUC.2UC.3) 

Also  store  all  and 


Kalman  filter  outputs. 


7 


■mm 


C.^  M  {i<l|i  s  ^  ♦)} 

=  set  of  previous  track  sets  which 
will  be  nontrivially  updated  by 
new  data,  given 
^  ♦?.  _ 


Go  to 


Cj4  t 


Set  C., 

3I 

and  C . n 
Go  to© 


-  2  log  pr(Z^^^^ 

(ieCj^UCj2UCj5UCj^)  . 


Also  store  all  L.'^-'-'*s  and 

1 

Kalman  filter  outputs 


Go  to  d) 


C- 


set  of  all  previous  track  sets  which  will 


be  updated  by  predictions  only  (no  new  data) , 

given  _ 

jL-42 - 


"Yes  I  No 

j__:  _ 1 


/Set  C.2  “ 
Yes  i 


^  Go  to  1^.  ^  C^4  /  jco  to 


No 


fSet  0-2  =  ^ 
and  c  ■  A  =  ^ 

34 

Go  to  (D 


_ ^ 

{ill<  i  §  $  (and  =  ,|,)} 

=  set  of  all  new  (one  point)  track  sets,  given 

q'^’  , _ _  —J 


Inputs;  For  i  e 


From  before:  X.  .  A.  s  s.  »  .Li""  "" 

A, 3  1, 3-1,3  1  *  A 

_ _  a=0.  . 


,  l.(M)“  L.O-l) 

'  1  9 


At  present:  r.j^  known  positive  number.  B..(r.,  by  m)  for  k  such  that  Z 

^  J  ^  J  ^  1  j 

*^3k  ^^jk  ^jk  positive  definite  ),  for  k  such  that  ^ 

^^ijkl^  ^  ^jl*  ^  ^  such  that  Z..j^  4,} 


®ii 


^  q  \ 

such  that] 

Ziik  J 


s/^5  ^  ,  0-1)  ,  3 

11  ij 


v-.-i  —  2.--1  “B.i'X..  .  ieC-- 

13k  13k  3k  13-1.3  j 


Go  to  ) 


^4  4=  y  B7,R:iB., 

13  ^  31^  jk  3k 

1  ^  k  ^  q 
such  that 

^ijk  ^  ^ 


A. . .  =  (D. .  +  A.  .  ,  .) 
133  13  i,3’1.3 


=  D.V  •  (I  -  (I  +  D.  .A.  .  ,  O'^) 
13  m  ''  m  13  1,3-1, 3''  '' 


*"  ^  1/’  ■^Ir  ^  ^4'L--'  *  ^4  1,  R^ 


T  ,.-1 


k'k^  ''jk'  ‘'jk"  3k'' 133  ^jk"  3k 


for  l^k",k"4q  such  that  ^  4),  4. 


log  det  E..  =  V  log  det  R., 
^  J  V  3  ^ 

(  such  that  ] 

V^ijk  *  *  J 

*  log  det  (I^  + 


Continue  track  set  i 
Redefine  Cjj  to  be 

restricted  to  i  such 

that  (g) , 

Yes  No _ _ 


M  is  a  fixed  large  number 
Begin  formally  new  track 
set:  replace  i  by  M+i. 
Put  M+i  e  C .  , 

set  Z(M’'  =  ♦ 

M+1 


Go  to  (D 


Go  to 


fGo  to 


C-12 


<X  I  I  <X  I  |<X 


© 


C-13 


m 


Inputs:  ieC. 


From  before:  = 

1 


Formally  set  at  present:  X.  =  E(X..),A.  .  ,  .  =  Covrx..) 

13  '•  xj-' 


At  present:  ^  such  that  f  ^  » 

by  T.^  positive  definite)  for  k 

such  that  Z-.,  ¥  ♦  , 
ijk  ^  ^  * 

{Zf^kl^  ^  ^j25  such  that 


I  such  that) 
Vijk  ^  ♦  y 


'’ijk  “  ^ijk’^jk^^^ij) 


For  each  i e  (for  given  new 

Let  4a  ^  j  ,  l^k^q  such  that 

^iak 

be  formed  out  of  at  least  two  old  track 

sets  or 

former  false  alarms.  Fix  j 

Let  be  the  minimal  value  of  a. 

for 

such  that  Z.^  ^  <j>. 

- - , - - - - — - - 

1 

In  ©.  replace  j 

everywhere  by  a.  Go  to 
and  loop  through  only  up  to, 
and  including, 


a+i<j  ?y 

^i,a+l  ^ 

1  No 


Go  to  0^ 
and  loop 
through,  up 
to,  and  in- 


Go  to 
and  loop 
through  up 
to.  and 


I  Outputs : 

7i,i  .  J+l» 

^  for  i  e  C  .  . 

1  j  4 

Return  to  j  notation 


i.e.,  no  previous  false  alarms 

are  destroyed  and  Zq.  ={Zq^j,  l<k4q}j<<J'? 


For  all  k,  1  ^  q 
such  that  ^  ♦: 

Inputs:  rjj^  known  number  ^  1 
known  r^j^  by  1 

M.,  known  positive  definite 

J  r.,  by  r., 
jk  ^  3k 

fjk  “  kno\%Ti 

'  ^Ojk  * 


^Ojk  "S  ^^Ojk«'®3k^^’^’jk^5o3ku*®jk^ 

U=1 


‘-Ojk  '  2  ”jk’ 


T  =  R  +  1  "' 

^Ojk  *^03  k  ^Ojk 


2  log  prCZQ^5)|Q(i)) 

=  -2  log  prCZg^j'^^lQ^j'^^) 


Z  (j)  =  7  (3-1)7 

i.e.,  =  t, 

but  Zq.  f  4,? 


♦  ^  Z, 


(3-1) 


Go  to  (35 


Go  to  (34 


-2  log  prCZg^i^ 

-2  log  pr(Zgf3-l)|Q(J-l) 


4-  2:  . 

1  ih  ;<  q\ 

such  that 

^Oik  ^  ^ 


^Ojk 


-2  log  pr(Z„^^^  |Q^^^,) 


E  Sjk 

1  ^k  ±q  \ 

such  that) 

'oik  *  ♦/ 


6 


everywhere  replace  j  by  a 


Is  o  such  that  there 
are  k's,  l^k<q 
such  that  #  ♦? 


No _ 

,  I _ 

-2  log  prCZp^®^ 1Q^“^) 

«  -2  log  pr(ZQ^“’^^|Qf“-^^ 


Replace  a  by  a+1 

=1:=; 

0+1 ^  j  ?  j 


Yes 


Loop  through 

(3^  up  to,  and 

including , 
- - 


Yes 


No 


Output : 

-2  log  prCZ^^^^ 


Return  to  3  notation 


Kon^Geolootlon  Target  Attribute  Data 

© 

(  Normal  Approximation  Used^ 


Yes 

r - - - - 

Go  to 

No 

(yCj-l! 

Yes 

No 

Thus,  «  Yj  ?<  ♦ 


For  all  i  such  that  Y 


21og  pr(Y^^^|Q^^^) 
2  log  pr(Y^^  Iq^ 


;  -^Wi:  IT  pr(Y,j, 

O  VH  <-r  /I  ^  V  ^  «  V 


Kalman  Filter 


For  all  i>l  such  that  e  C  , 


pr(Y,0)|„.,Qa))  . 

I  such  that) 

Vi«k  ’■  ♦  / 


=  -2  log  2  pr(Y.  -prCH.) 

1  u  ^  ■‘■ 


H^cC 


-  2  log  pr(Y^^^  = 


E 

l^i  \ 

such  that! 

Y^Cj)  y 


Kalmen  Filter 


Approach  for  Non-geolocation  Target 


Attribute  Data  Under  the  Normal  Distributional  Approx imat ions 


Go  to  and  everywhere  replace: 


1  by  b 

’"jk “ijk 

■ijk  ’^ijk 

"^jk  '*ijk 

by  yO) 

®jk  ®ljk 

by 

♦j  '>>'  Ib 

etc. 

G.  by  0 

J 

.  by  L 

by  0 

by 

{  Thus  27  becomes 

■ij  ‘■ij 

H. .  .  T  =  H. . . 

iJJ 

“  A***  ) 

11  »j+i  111  * ' 

D. .  by  V. . 

11  '  13 

etc. 

. .  by  H. , 

13  ^  1  ’ 

^ijk'  ,k"  ^ijk'k" 

A 

.  .  ,  .  by  H.  .  ,  . 
i;3-l0  ^ 

mm 

Z. .  by  Z . . 

11  ^  11 

A 

.  .  .  by  H.  .  . 

'  I'.IO 

Sjk'k**’’'  ®ijk'k' 

etc. 

F. .  by  F. . 

11  ^  13 

etc. 

Also  everywhere  in  replace: 

by 

"ij  "IJ 

^ _ njii  ''ijk- _ 

Further  evaluations: 

For  use  (4.98). 

If  >b: 

For  Sjj .  use  eqs.  (4.93)  -  (4.95). 

For  ,  use  (4.91),  (4.101),  (4.102). 

For  log  det  >  use  (4.103). 

Si^jk  *^ijV^ijk  *  (4.96),  (4.97). 

For  <b: 

^  'J' 

For  Sj^jk"  »  analogue  of  eq.  (4.102). 

®/jk'*^ijk"lr8ijlr  *  ose  eqs.  (4.109)  -  (4.111). 

®ijk^’''(ijk"k-'‘*ijk-'  »  ose  eqs.  (4.98),  (4.109), 

(4.112),  (4.113). 

When  is  reached,  skip  to  and  continue  looping  through 

as  usual. 

When  is  reached  and  outputs  deposited,  resume  old  notation 
again  in  until  next  cycle  into 


Direct  Regression  Approach  to  Non-geolocation  Target  Attribute 


Data  ,  Under  Normal  Distribution^  Approximations 


log  det 


-  Z  H. 


1  j 


5^  log  det  f  R  ) 


log  det  =  2^  log  det  R^^j^ 

^  /^<k<q\ 

/  such  that] 


log  det  R 


a . 

iak  “ 


log  r, 

°iok7 


'i  di  '  {Vr)  •  £  <  1  ) 

'  V'  /li,T<q,  3l  <«<j. 


such  that  ^  4>, 

and  there  is  n,  l^n^b 
such  that  Y=b.^j^^ 

^  ^iak  y 


log  det  (I^  +  P^^^-Cov(H^)) 


«  log(l  +  d/^^^-varCH^^)) 

Y  =  1 


log  det  Cov(Y^^^^ 

«  log  det  R^^^  +  log  det  (Ij^  +  ^ 'CovCH^)) 


log  2it  +  log  det  Cov(Y^^j^ 


c^Cj) 

®i»Y 


(Vr  )• :  2  £  (Vi„^  - 

/  such  that  ^  \ 

I  and  such  that  there  J 
\is  n,  l^ri^b,  so  that/ 


‘i.b, 


iakn 


^"^iakn  '^iotk 


Approximation  desired  minimizing 
knowledge  of  E(H^),  Cov(H.)  ? 


Yes  No 
lutation 
below. 


such  that 


Y.  1, 

s.  lak 


m 


SUMMARY  OF  FLOW  CHART  INPUTS  OUTPUTS 
Summary  of  T ime - Independent  Quantitites 

C  is  a  known  set  of  possible  non-geolocation  target  attributes. 
Card  (C)  =  a.  Each  element  of  C  is  a  b  by  1  vector  of  values 
such  as  (true)  hull  length,  flag  color,  identification,  etc. 

Fixed  number  of  sensor  systems  operating  is  q,  though  some 
at  times  may  fail  to  collect  data. 

m  is  the  common  dimension  of  each  target  state  parameter  vec¬ 
tor;  n  is  the  dimension  of  the  driving  noise  vector  Ceq.  (2.1)), 

If  normal  approximations  are  used  for  the  non-geolocation  tar¬ 
get  attributes: 

E(H^),  Cov(H^)  are  the  mean  and  covariance  matrix  of  variability 
or  randomness  for  any  target  i  ,  and  it  is  assumed  their  values  are 
not  dependent  on  i  (nor  on  j ) . 

If  normal  approximations  are  not  used  for  the  non-geolocation 
target  attributes: 

pr(H^)  the  prior  probability  function  of  e  C  is  known; 
it  is  usually  assumed  to  be  uniform:  pr(H^)  =  Va  for  all  e  C, 
where  a  *  card(C).  (pr(H^)  does  not  depend  on  j . ) 

Summary  of  Previous  Outputs  for  Data 

Partitioning  •••}»  where 

is  the  perceived  (target,  if  i^l,  false  alarm,  if 
i  =  0)  i^^  track  set  consisting  of  data  *  {Z.  v  I  0 < a  « j , 

l<k<q,  such  that  2.  ^ 

“  “  lot  K 
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z> 


Kalnan  filter  outputs:  X..  ,  A.  .  ,  . 

(m  by  1)  (hi  by  m  positive  definite) 


(b  by  1)  (b  by  b  positive  definite) 


.  (The  last  two  are  not  re 


quired  if  the  Kalman  filter  -  normal  approximation  procedure  is  not 

used  for  the  non-geolocation  target  attributes.) 

( 

=  -2  log  pr(Z^^j"^^  for  i>0;  ,  for  i>  1; 

=  -2  log  pr(Y.  IQ^^'^^),  for  i>l; 

-2  log  pr(z  =  S  TL.tM))  ; 

i=l,2,..  ^ 

-2  log  pr(Y  )  ; 

=  -2  log  pr(Z^^j'^^ 

-2  log  pr(ZQ^j’'^^lQ^^'^)) 

-2  log  pr(Y^3‘^^ 


Summary  of  Data  Inputs  at 
Sampling  Time  t ^  but  Before  New 
Partitioning  Q^^^is  carried  out 

Consider  sensor  system  k,  k  =  1,..,  q: 

If  no  new  data  at  all  arrives,  Zj  =  4» 

If  some  new  data  arrives  Z.  7*  and  for  each  sensor  k 

•  J 

which  obtains  data,  m^j^  ^1  is  the  number  of  data  reports  received; 
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^3  ■  ^*V.3.k’ 


1 ,  .  .  jilij 

k=l y . . 


(where  above,  set  =  41  for  those  k*s  for  which  =  0). 

The  report  received  by  sensor  k  at  t^  is  , 


2yjlf  is  rjjj  by  1  geolocation  data  vector 


t-'jk  i  1’ 


vector 


Z'jj^  is  a^jj^  by  1  non-geolocation  data  target  attribute 


=  ^Yjk  »  ^  ' 


It  is  possible  for  a  given  y,  l^Y^^j^  ^  for  1,  to  be 

such  that  ^  OY  ~  both  (otherwise 

”jk  • 

Zy-y.  -  ♦  iff 


Yjk 


0. 


Z'jk  contains  as  components,  typically  measured  positions, 
velocities,  etc.,  from  sensor  system  k  at  tj  for  report  Y 

Zyjk  Contains  as  components,  a  given  subset  (of 

size  of  a  fixed  set  of  b  attribute  values. 


C-27 


iJP,  IJIiL;mWIIJ|U. 


Typically,  this  could  be 


'Yjk 


(number  of  radars  on-board\ 
hull  length  I 

flag  color  / 


where  the  total  attribute  set  is  {number  of  radars  on-board, 
identification,  shape,  hull  length,  flag  color} 

{b  «  5  here) , 

Also  given  for  eapJr  Sensor  k  is  ®jk  ^  ’'‘jk  measurement 

matrix. 

Typically,  if  the  i^^  target  state  parameter  vector  is 


(m  »  4) 


9 


and  only  x-  and  y-  positive  measurements  can  be  made  at 
by  sensor  k,  ...then 


(observed  x-pos.  at  t ^ \ 
observed  y-pos.  at  t./ 


and 


®ik 


^1  0  0  o\ 

V 0  1  0  0/ 


(2  by  4) 


(r^k  -  2)  , 


(see  eq.  (2.2)) 


Also  given  is  the  measurement  error  covariance  matrix. 

*^jk  (tjk  positive  definite)  is  typically  obtained 
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T 

approximately  as  *  where  (f is  a  matrix 

of  Jacobians  of  the  transformation  from  sensor  coordinate  to 
cartesian  coordinate  space,  and  is  the  covariance  matrix  of 

measurement  error  of  sensor  system  k  at  tj  for  its  natural 
(sensor)  coordinate  system.  (See  also  eq.  (2.2).) 

Also,  for  each  seftsor  system  k,  is  given  (^jk  ’ 

the  mean  of  the  false  alarm  dispersion  ;  the  corresponding 

(rjk  ^jk)  positive  definite  covariance  matrix. 

Also,  the  common  target  motion  transition  matrix  fj 
(m  by  m)  driving  noise  coefficient  matrix  Gj  (m  by  n)  and 
driving  noise  covariance  matrix  Pj  (n  by  n  positive  definite) 
are  all' known.  (See  eq.  (2.1).) 

If  required  (usually  by  a  fixed  procedure  ),  E(X^j)  (m  by  1) 
and  Cov(X. .)  (m  by  m  positive  definite)  will  be  known,  but  not 
dependent  on  i . 

Similarly,  for  non-geolocation  attribute  data,  given  are: 


B 


Yjkl 


B 


Yjk 


B 


Yjka 


Yjk 


®Yjkn  “  CO,. . ,  0,  1,  0,  . . ,  0) 


where  the  1  occurs  in  the  b 


th 


Yjkn 


position,  1  4  n  4  ^Yjk’ 


where  1  <  i  b 


and 
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's  are  subvectors 


0 


Quantities  Known  at  t j 


Only  After  Partitionint 


is  Carried  Out 


•••  y  known,  where 

fjjj  the  number  of  false  alarms  for  sensor  system  k  at  tj 

is  known  and  index  i  ’representing  target  i  is  now  known;  all 
relative  to 

Following  rearrangements. 


i=0 ,1,2,... 
k*i, . . ,q 
for  t  ♦ 


Vijk/  VSjk/ 


for  some  y  »  Y  =™jk* 

Zfjk  is  r^j^  by  1  geolocation  data  vector  assigned  to 
track  set  i 

^ijk  ^ijk  ^  geolocation  data  target  attribute  vector 

assigned  to  track  i 

^ijk  ^ijk  vacuous;  ^  always. 

Thus,  become  determined;  similarly  for 

prCYi-^lHf). 


The  total  data  Z up  to  t^  can  be  broken  up  two  different 


ways : 
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ilMLli  ^  ,J|  u. 


where 


-  iZ^^K  V0)> 

f  * 


-  0,  1,  2.  ...} 


(j)  -  {Z.^j^li  »  1,  2,  ...}  , 


2i 


0) 


(Zi^l^fo  <  a  <  j;  l<k<q. 


geolocation  data  for  track  set  i  up  to  tj  , 


(3) 


^iok  ^ 


false  alarm  data  up  to  t^  , 

-  1,  2,  ..)  , 


’fl 


(3)  . 


<^iakl  »  i  “  S  3  ;  1  <  k  <  q,  +  ♦) 

■  non-geolocation  attribute  data  for  track  set  i 
up  to  tj 

In  turn,  these  determine  sets  Cj2»  Cj^,  C-j,  Cj2 

C^j  Cj^  (the  last  four  are  not  required,  if  the  Kalman  filter- 
normal  approximation  procedure  is  not  used  for  the  non-geolocation 
target  attributes).  (See  Flow  boxes  7-10  for  definitions.) 
Note  also  the  notation: 


f.  . 

13 


^^ijk^ 


such  that 
^ijk  ’*  ♦ 


;  similarly  for  Y 


13 
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